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if 


Enclosed  please  find  the  1986  Gravity  Gradiometer  Conference  Presentations. 

I  would  like  to  apologize  for  taking  my  time  in  getting  this  out  to  you. 

I  know  a  year  late  is  a  long  time  after  my  promised  April  86  delivery. 

For  those  of  you  who  are  attending  this  years'  conference,  please  remind 
me  of  how  long  it  took,  so  I'll  be  embarrassed  enough  to  get  it  out  more 
promptly. 

The  experiment  I  tried  with  the  “Question/Answer  Sheets”  was  a  partial 
success.  However,  as  with  any  experiment  there  was  experimental  error. 

Some  data  points  are  missing,  which  leads  me  to  believe  that  I  did  not 
stress  the  importance  of  experiment  sufficiently  to  everyone.  To  all 
those  speakers  who  may  feel  that  their  responses  were  not  accurately 
recorded,  you'll  have  to  blame  me.  It  is  my  hope,  though,  that  no  bias 
was  introduced  into  the  experiment.  It  was  my  conclusion  that  most 
errors  were  of  a  random  nature  and  have  not  detracted  from  the  resulting 
compilation.  You  will  also  have  to  blame  me  for  typo's,  spelling, 
grammatical  errors  and  any  other  mistakes  that  went  uncorrected.  In  some 
cases,  I  have  altered  the  wording  of  the  questions  or  responses  for 
clarity. 

Thank  ^ou  for  yo/r  participation  at  the  14th  Gravity  Gradiometer  Conference. 

i'lJ-Lt 

— J^^dalt,  Capt,  USAF 
Geodesy  and  Gravity  Branch 
Earth  Sciences  Division 
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Fourteenth  Gravity  Gradloaetry  Conference 
United  States  Air  Force  Academy 
Colorado  Springs,  Colorado 


February 

-  Buses  depart  USAFA  Officers'  Club  for  Fairchild  Hall 

-  Registration  -  3rd  floor  Fairchild  Hall,  South  End 

-  Welcome/Introduction  -  ILt  Terry  J  Fundak 

-  Opening  Remarks  -  Dr.  Donald  H.  Eckhardt 

-  "  Development  Experience  of  a  Moving  Base  Gravity  Gradiometer  and 

Discussion  of  Future  Applications*' 

Mr.  Ernest  H.  Metzger 

-  "  Gradlometry  &  Geodesy,  or  Separating  Inseparables  " 

Dr.  Christopher  Jekell 

-  "  Requirements  for  the  use  of  Airborne  Gradlometry  in  Exploration 

Geophysics  " 

Dr.  Klaus-Peter  Schwarz 

-  Break 

-  "  Applications  of  Superconducting  Gravity  Gradiometer  System  toward 

Inertial  Guidance  and  Fundamental  Science  " 

Dr.  Hlnghung  A.  Chan,  Dr.  Martin  Vol  Moody,  Dr.  Ho  Jung  Paik* 

-  "  Quick  Review  of  Gradlometer-Alded  Land  Navigation" 

Dr.  Warren  G.  Heller 

-  "  Efficient  Gravity  Gradient  Data  Gathering  " 

Mr.  Michel  Bllello*,  Dr.  John  B.  Breakwell,  Dr.  Daniel  B.  DeBra 

-  "  Use  of  Terrain  Elevation  Data  in  Airborne  Gradlometry 

Dr.  Julian  L.  Center,  Jr. 

-  Buses  depart  Fairchild  Hall  for  USAFA  Officers'  Club 

-  Luncheon  -  Officers'  Club 

-  Buses  depart  Officers'  Club  for  Fairchild  Hall 

-  "  The  Gravity  Gradiometer  Survey  System  " 

ILt  Terry  J.  Fundak 

-  "  Airborne  Gravity  Gradiometer  Data  Processing 

Dr.  William  J.  Hutcheson 


r 


1410  -  ”  Isomorphic  Geodetic  and  Electrical  Networks:  An  Application  to 
the  Analysis  of  Airborne  Gravity  Gradloneter  Survey  Data  ” 

Dr.  Donald  H.  Eckhardt 

1430  -  "  Stage  II  Processing  of  Airborne  Gravity  Gradlometer  Data  using 
Frequency  Domain  Techniques  ” 

Mr.  Anthony  A.  Vasslllou 

1450  -  Break 

1505  -  "  Karhunen-Loeve  Gravity  Gradlometer  Data  Processing  " 

Dr.  Sam  C.  Bose 

1525  -  "  Gravity  Gradlometer  (GGSS)  Test  Planning  and  Test  Data  Treatment 
Dr.  Warren  G.  Heller 

1545  -  "  Gradient  Integration  Procedure  for  Path  Error  Reduction  " 

Mr.  Alan  E.  Rufty 

1615  -  Buses  depart  Fairchild  Hall  for  USAFA  Officers'  Club 

1630  -  Reception  -  USAFA  Officers’  Club 


♦ 


Wednesday  12  February 


0800  -  Buses  depart  USAFA  Officers'  Club  for  Fairchild  Hall 


0810  -  "  NASA  Requirements  for  a  Spaceborne  Gravity  Gradlometers  -  An  Overview 
Mr.  Charles  J.  Finley,  Dr.  David  E.  Smith 
Presented  by  Mr.  Werner  D.  Kahn 


0825  -  ”  Gravity  Field  Fine  Structure  Mapping  using  a  Spaceborne 
Gravity  Gradlometer 
Mr.  Werner  D.  Kahn 


0845  -  "  Superconducting  Gravity  Gradlometer  on  the  Space  Shuttle  " 

Dr.  Samuel  H.  Morgan,  Mr.  Joe  R.  Parker* 

0910  -  "  Development  of  Superconducting  Gravity  Gradlometer  for  Space 
Applications  " 

Dr.  Hlnghung  A.  Chan,  Dr.  Martin  Vol  Moody*,  Dr.  Ho  Jung  Paik 
0930  -  Break 

0945  -  "  Platform  Requirements  and  Error  Compensation  for  a  Superconducting 
Gravity  Gradlometer  ” 

Dr.  Ho  Jung  Pa Ik 

1015  -  ”  TLC  for  a  Magnetically  Floated  Gravity  Gradlometer  " 

Dr.  Dave  Sonnabend 


1040 


Development  of  a  High-Sensitlvity 
Gradlometer  for  Spaceborne  Use  ^ 
Dr.  F.  Bordonl,  Dr.  F.  Fullglni*, 
Dr.  Enrico  C.  Lorenzinl 


Non-Cryogenlc  Gravity 


Dr.  B.  V.  lafolla. 


1110  -  “  Coaoon  Mode  Balancing  Gradioaeter  with  Monocrystalline  Silicon 
Suspension  for  High  Sensitivity  Gravity  Measurements  ” 

Dr.  Jean-Paul  Richard 

1140  -  "  Liquid-Supported  Torsion  Balance  as  a  Gradiometer" 

Dr.  James  E.  Faller,  Mr.  Paul  T.  Keyser* 

1210  -  Closing  Remarks 

ILt  Terry  J  Fundak 

1230  -  Buses  leave  for  lunch  at  USAFA  Officers'  Club 
1345  -  Tour  of  Air  Force  Academy 
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DEVELOPMENT  EXPERIENCE  Cf  A  MOVING  BASE  GRAVITY  GRADIOMETER 
AND  DISCUSSION  OF  FUTURE  APPLICATIONS 

E.  H.  Metzger 
Bell  Aerospace  Textron 
P.  0.  Box  One 
Buffalo,  NY  114240 


ABSTRACT 

A  sunmary  of  the  development  experience  of  the  rotating  accelerometer 
gravity  gradiometer  from  its  conception  to  reduction  to  practice  yielding 
accurate  measurements  aboard  a  moving  vehicle  is  presented.  Potential  future 
applications  of  moving  base  gravity  gradiometers  are  outlined  and  the  techni¬ 
cal  difficulties  to  achieve  these  objectives  are  discussed. 
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SCHEMATIC  ILLUSTRATION  OF  ROTATING  FIXTURE 


TYPICAL  ACCELEROMETER 


TIME- HOUR 


QGI  SIGNAL  PROCESS  BLOCK  DIAGRAM 


MASS  ATTRACTION  TEST  -  MODEL  VM  ACCELEROMETERS.  FIBERGLASS  COATED  FLASK  BETWEEN  ACCELEROMETERS 


a> 


RECORDING  OF  DATA  OBTAINED  DURING  MASS  DETECTION  EXPERIMENT 


RESULTS  OF  MASS  DETECTION  EXPERIMENT 


TEXTRON 


COMPARISON  OF  THEORETICAL  AND  MEASURED  THERMAL  NOISE 
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ROTATING  ACCELEROMENTER  GRAVITY  GRADIOMETER 
ON  HOLLOMAN  RATE  TABLE 


MOUNTING  BLOCK  DECOUPLING  INTERFACE 
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INSTRUMENT  BLOCK  AND  ELECTRONICS 
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LONG  WAVE  LENGTH  GRAVITY  DISTURBANCE  VECTOR  DETERMINATION 


GRAVITY  GRADIENT  SIGNAL  <E 


HARMONIC  DEGREE  (S) 


Bell  Aerospace  CS3ZIIEID 


r<* 

CM 


CO 

Ui 

O 

UI 

CO  H 

M 

t  g 

5  ? 

CO  ^ 


=> 
5  O 

^  I 
CO  ^ 

<  s 


u. 

o 

i  » 

^  i 

UI  CO 

E  5 


1*5 

UI  o 

S  ^ 

5  u. 
=  2 
U  O 

>  P 

CD  O 

>  ^ 
j  a 
u.  u- 

UI  < 

t.  O 


<  b 

CO  Q 

£  8s 

H  Z 

UJ 

«  UJ 

is 

1“ 

^  h- 

H  UJ 

u 

<  H 


U  UI 

u  u 

CO  Z 

O  S 

“  o 

CM  U. 

UJ  tc 
*>  UJ 

o  u- 
d  M 
u.  X 
O  I- 

-J  UI 
UI  O 

^  i 

-I  o 

s  % 

i  § 

K  O 

UI  O 
H  flc 

tit  V 


O  UJ 

<  2 
a.  Q 
CD  < 
<  oc 

If 

*-  < 
UJ  ^ 
HI  tti 

u!  o 

IL  < 

O  u* 

z  2 

UJ  CD 


V 

00 


SATELLITE  OR  REENTRY  BODY  VERTICAL  ERECTION  SYSTEM 


8 


oc 

o 

> 

o 

s 

> 

oc 

H 

Z 

UJ 

Ul 

oc 

< 


H  OC 
OC  < 

u. 

>o 

LLl 

P 

ii 

o< 
^  o 

D  H* 

8  “ 

w  UJ 

OC  z 

1“ 

OC  UJ 

a  M 

>  UJ 

H  »- 

>  n 

nr  Ul 


s 


n 


z 

UJ 

5 

< 

oc 

o 


> 

H 

> 

< 

OC 

o 


< 

z 

i 

o 


UJ 

z 

H 


O 

S 

V) 

UJ 

O 

3 

H 


< 

< 

H 

s 

oc 

o 


d 

UJ 

OC 

5 

o 

UJ 

OC 

UJ 

00 

> 

< 


00 

oc 

o 

z 


oc 

UJ 

K 

Ul 

z 

o 

5 

< 

oc 

o 

UJ 

z 

H 


O 

z 

o 

p 

< 

C 

00 


< 

u 


CO 

< 

S 


Bell  Aerospace 


OIL  EXPLORATION  THROUGH  BORE  HOLE  GRAVITY  GRADIENT  MEASUREMENTS 
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GSS  WILL  GIVE  STRONG  GRADIENT  RATE  OF  CHANGE  SIGNALS  TO  APPROACHING 


TERRAIN  AVOIDANCE 


PAPER  TITLE: 


DEVELOPMENT  EXPERIENCE  OF  A  MOVING  BASE  GRAVITY  GRADIOMETEW 
AND  DISCUSSION  OF  FUTURE  APPLICATIONS 


Speakers  name: 


Ernest  H.  Metzger 


Questions  and  Coniments  : 

Milton  Trageser:  I  have  never  seen  data  which  show  your  Gradioraeter 's  quick  time 
response  and  noise.  Can  you  show  some  data  on  this? 

SPEAKERS  RESPONSE:  (No  Response  Noted) 

Daniel  DeBra:  What  is/are  the  principle  hardware  differences  in  the  Navigation 
GGS  and  the  Mapping  GGS? 

SPEAKERS  RESPONSE:  The  Mapping  GGS  is  more  sensitive.  The  accelerometers 
have  fewer  turns  on  the  proof  mass  to  require  more  A/g  (Amoeres  per  unit  of  cravitv). 
This  has  produced  a  significantly  better  signal  to  noise  ratio.  It  is  not  needed  for 
the  navigation  application. 


0 


GRADIOMETRY  AND  GEODESY,  OR  SEPARATING  INSEPARABLES 


C.  Jekell 

Air  Force  Geophysics  Laboratory 
Earth  Sciences  Division 
Hanscom  AFB,  MA  01731-5000 


ABSTRACT 

The  purpose  of  geodesy  is  to  determine  the  size  and  shape  of  the  earth  and  to  deter¬ 
mine  the  exterior  gravity  field.  These  two  facets  of  geodesy  are  tightly  interwoven. 
In  fact,  a  knowledge  of  the  gravity  field  is  a  prerequisite  for  terrestrial  position¬ 
ing,  and  vice  versa,  in  the  most  comprehensive  solutions  (Integrated  geodesy),  both 
j^re  determined  simultaneously  from  all  available  measurements. 

Two  methodologies  being  proposed  for  rapid  local  gravity  mapping  are  airborne  gravime¬ 
try  and  airborne  gravity  gradlometry.  In  addition  to  the  classic  chicken-egg  problem 
mentioned  above,  there  is  the  problem  that  airborne  gravity  instruments  sense  not 
only  gravity  (gravity  gradients)  but  also  various  effects  of  an  accelerated  coordinate 
frame. 

A  rudimentary  analysis  shows  how  gravity  gradlometry  mitigates  these  problems  of 
inseparability. 


GRADIOMETRY  8  GEODESY,  SEPARATING  INSEPARABLES 


PROBLEMS  OF  MOVING-BASE  GRAVITY  MAPPING: 

-  INSTRUMENTS  IN  AN  ACCELERATED  FRAME 

-  GRAVITY  FIELD  IS  A  FUNCTION  OF  POSITION 

SIMPLIFIED  ANALYSIS: 

-  USE  SPHERICAL  APPROXIMATION 

-  LINEARIZE  OBSERVATIONAL  MODEL 

-  ASSUME  LOCAL  LEVEL  STABILIZATION 

-  ASSUME  NO  HORIZONTAL  ACCELERATIONS,  CONSTANT  ALTITUDE 

-  GRAVITY  (GRAVITY  GRADIENT)  DETERMINATION  AT  ALTITUDE  TO 


ACCURACY  OF  1  MGAL  (1  E) 
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REQUIREMENTS  FOR  THE  USE  OF  AIRBORNE 
GRAVITY  GRADIOMETRY  IN  GEOPHYSICAL  EXPLORATION 

K.  P.  Schwarz 
A.  A.  Vasslllou 
The  University  of  Calgary 
Division  of  Surveying  Engineering 
2500  University  Drive  N.W. 

Calgary,  Alberta  T2N  1N4 
CANADA 


ABSTRACT 

Gravity  methods  In  geophysical  exploration  exploit  the  fact  that  density 
^rlatlons  In  rocks  give  rise  to  variations  In  the  anomalous  gravity  potential. 
Interpretation  of  these  variations  gives  Information  on  rock  formations  that 
may  be  of  geophysical  Interest.  The  anomalous  gravity  potential  of  some 
typical  formations  are  specified  In  terms  of  wavelengths  and  amplitude  and 
the  requirements  for  resolving  these  signals  by  single  axis  or  multi-axis 
gravity  gradlometers  are  stated.  Simple  numerical  examples  are  used  to  show 
the  effect  of  downward  continuation,  aliasing,  and  profile  spacing.  Finally, 
the  noise  spectra  of  current  gravity  gradlometers  are  discussed  with  view  to 
their  potential  use  In  these  applications. 


REQUIREMENTS 
FOR  THE  USE  OF  AIRBORNE 
GRAVITY  GRADIOMETRY  IN  GEOPHSICAL 

EXPLORATION 

by 

•  K.P.  SCHWARZ  AND  A.A.  VASSILIOU 

UNIVERSITY  OF  CALGARY 


GRAVITY  METHODS  IN  EXPLORATION 


Major  exploration  areas 

•  PETROLEUM  EXPLORATION 

•  MINERAL  EXPLORATION 

•  WATER  RESOURCE  EXPLORATION 

Major  exploration  methods 

•  SEISMIC  (petroleum) 

•  ELECTROMAGNETIC  (MINERAL) 

•  MAGNETIC 

•  GRAVIMETRIC 

Gravity  methods 

•  COMPLEMENT  SEISMIC  AND  ELECTROMAGNETIC  METHODS 

•  ARE  INEXPENSIVE 

•  GIVE  MULTIPLE  SOLUTIONS  (NON-UNIQUENESS  OF 
POTENTIAL  methods) 

•  GIVE  LITTLE  GEOMETRIC  INFORMATION 

•  ARE  SLOW 

Advantages  of  airborne  gravity  gradiometry 

•  MORE  COMPLETE  INFORMATION  (FULL  TENSOR) 

•  HIGHER  RESOLUTION  POSSIBLE 

•  FASTER  DATA  ACQUISITION 

•  LARGER  AREAS 

•  IDEAL  COMBINATION  WITH  AEROMAGNETICS 


Integrated  approach 


APPROACH 


Requirements  in  terms  of  wavelength  and  amplitude 

•  TYPICAL  STRUCTURES 

•  RESULTING  WAVELENGTHS  AND  AMPLITUDES 

Signal  detection  at  flying  altitude 

•  attenuation  of  second-order  GRADIENTS 
WITH  ALTITUDE 

•  DETECTABILITY  OF  OFF-PROFILE  SOURCES 

•  GEOLOGICAL  AND  TOPOGRAPHICAL  NOISE 

•  SAMPLING  RATE  AND  MEASURING  ACCURACY 

Downward  continuation 

•  AMPLIFICATION  OF  NOISE 

•  RECOVERY  OF  DESIRED  SIGNALS 

Use  OF  CURRENT  AND  FUTURE  SYSTEMS 

•  CURRENT  ROOM  TEMPERATURE  SYSTEMS 

•  FUTURE  ROOM  TEMPERATURE  SYSTEMS 

•  CURRENT  SUPERCONDUCTING  SYSTEM 


FUTURE  SUPERCONDUCTING  SYSTEM 


OIL  TRAPS  IN  SEDIMENTARY  STRUCTURES 


SALT  DOME  LIMESTONE  REEF 


GAS 


FAULT  TRAP 


UNCONFORMITY  TRAP 


TYPICAL  WAVELENGTHS  AND  AMPLITUDES 


Structure 

Range 

•  • 

Wavelength  Amplitude 

Salt  dome 

3  -  10  KM  5  -  30  E 

Reef 

0.5-  5  KM  0.3-  5  E 

Anticline 

3  -  15  KM  3  -  10  E 

Fault 

1  -  5  KM  1  -  30  E 

PiNCHOUT 

1  -  5  KM  1  -  5  E 

Unconformity 

1-5km  1-lOE 

Ore  body 

0.3 -0.5  KM  0.5-100  E 

Gravity  Methods 
Useful? 

Yes 

Often 

Yes 

Yes 

Only  if  density 

CONTRAST  IS 
SUFFICIENT 

Yes 


*  Signal  amplitude  and  wavelength  at  the  surface 


OF  the  earth. 


STRONG  ATTENUATION  WITH  ALTITUDE 
Example  :  Ore  body 


Depth  i  200  m  below  ground 
r  a  d i u  9 -  1 OOm 

density  contrast  =  1.0  gr/cin**3 
flying  altitude  -  0.0  m 


Depth  *  200  (B  below  ground 
r  a  d I u  s= 1 OOm 

density  contrast  =  1.0  gr/cm**3 
flying  altitude  =  200.0  m 
Profile  *  across  center 


NEAK  ATTENUATION  MITH  ALTITUDE 
Example  :  Fault  trap 

Example  <  Fault  trap 
tength^S^O  km  t  depth  ^  1*0  km 
density  contrast  *  0-2  gr/cm**3 
flying  altitude  ^  0*0  km 


-20  '10  0  10  20  30(ltm) 


Example  *  Fault  trap 
length«5.0  km  »  depth  *  1.0  km 
density  contrast  •  0.2  gr/cm**3 
flying  altitude  =  0.2  km 


-20  -10 


POOR  DETECTABILITY  OF  OFF-PROFILE  SOURCES 
E)(ampue  :  Ore  body 

Depth  *  200  m  below  ground 
r  a  d I u  5= 1 00m 

density  contrast  ~  1.0  Qr/cm**3 
flying  altitude  =  200*0  m 
Profile  *  250  m  of  the  center 


"6  -4  -2  0  2  4  6  (1<m) 

Depth  t  200  m  below  ground 
r  a  d I u  s  =  1 00m 

density  contrast  *  1.0  gr/cm**3 
flying  altitude  -  200*0  m 
Profile  *  500  m  of  the  center 


-6  -4  -2  0  2  4  6  (km) 

Note:  Larger  structures  are  much  better  to  detect.  Sources  at  a  horizontal 

DISTANCE  OF  FIVE  TIMES  THE  SOURCE  DEPTH  CAN  STILL  BE  SENSED, 


DOWNWARD  CCNTINUATION 


* 

Figure  gives  the  amplification  factor  for  wavelengths 
RANGING  FROM  ABOUT  10  KM  (1)  TO  0.^  KM  (20) 


downward  conilnuaiion  operator  of  Tzz 
flying  altitude  =  200.0  m 
grid  spacing  =  300.0  m 


20  -j 
16  ^ 
12 
8  ~ 
4  - 


|n  '(  I  I  r  T 1  I  I  I  I  I  I  I  I  I  I  I  I  I 


0.0  0.4  0.81.21.6  2.0  2.4 


(cycles/km.) 
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SIGNAL  DETECTION  AT  FLYING  ALTITUDE 


Attenuation;  Strong  attenuation  of  short  wavelength  signals 
(  <  0.5  km)  with  altitude,  If  a  specific  frequency  band  is 

WANTED^  FLYING  ALTITUDE  CAN  BE  USED  AS  A  LOW-PASS  FILTER. 

Detection  of  off-profile  signals;  Mixed  gradients  are  useful 

TO  DETECT  SIGNAL  SOURCES  WHICH  ARE  NOT  DIRECTLY  BELOW  THE 
FLIGHT  PROFILE.  In  GENERAL^  SIGNALS  CAN  STILL  BE  DETECTED  AT 
A  HORIZONTAL  DISTANCE  THAT  IS  TWICE  THE  DEPTH  OF  THE  SOURCE. 

This  would  favour  a  full  tensor  gradiometer  over  a  one  axis 

GRADIOMETER  (122)* 


NOISE  REDUCTION 

Geological  noise;  Cause  by  'unwanted'  density  changes  in 

the  upper  CRUST.  WAVELENGTHS  AND  AMPLITUDES  ARE  SIMILAR. 

Elimination  only  possible  by  joint  interpretation  of  all 

RELEVANT  DATA. 

Topographical  noise;  Caused  by  surface  topography.  Airborne 

MEASUREMENTS  MUST  BE  REDUCED  USING  DEM  AND  DENSITY  ASSUMPTIONS. 

Smaller  features  can  be  filtered  out  by  flying  altitude. 

Residual  systematic  effects  have  to  be  expected  but  will  usually 

BE  SMALLER  THAN  GEOLOGICAL  NOISE. 


SAMPLING  RATE  AND  MEASURING  ACCURACY 


Assumptions  :  10  data  points  per  smallest  period 

Flight  speed  :  100  knots 
Flying  height  :  200  m  (660  ft) 
Signal  to  noise  ratio  ;  5:1 


Structure 

Sampling  Rate 
(m)  (s) 

Measuring 
Accuracy*  (E) 

Salt  domes 
Anticlines 

300  - 1500  6  -  30 

4  -  7 

Ore  BODIES 
Pinnacle  reefs 

30  -  50  0.6  -  1.0 

o 

1 

o 

All  others 

100  -  500  2  -10 

-i 

1. 0-2.0 

*  TO  resolve  minimum  amplitude  at  flight  level. 


DOWNWARD  CONTINUATION 

Strong  amplification  of  noise  and  short  wavelength  features 
MAKE  the  interpretation  OF  THE  DOWNWARD  CONTINUED  SIGNAL  MUCH 
MORE  DIFFICULT,  SIGNAL  DETECTION  AND  INTERPRETATION  AT  FLYING 
ALTITUDE  ARE  THEREFORE  PREFERABLE.  AVAILABLE  GROUND  GRAVITY 
CAN  BE  UPWARD  CONTINUED  TO  FACILITATE  LONG  WAVELENGTH  DETECTION. 


USE  OF  CURRENT  AND  FUTURE  SYSTEMS 


Current  room  temperature  (Bell  system) 

Noise:  9  E  hz”^  ^  ,1  Hz  averaging 

- >  Large  salt  domes  and  anticlines 

Flying  altitude  of  600  m  is  acceptable  for  these  features. 

Future  room  temperature 

- >  lower  noise  level  with  room  temperature  systems  seems 

TO  BE  FEASIBLE.  HOWEVER^  NOISE  REDUCTION  HAS  TO  BE 
ACCOMPANIED  BY  A  HIGHER  SAMPLING  RATE>  TO  MAKE  THEM 
ADAPTABLE  TO  THESE  APPLICATIONS. 

Current  superconducting  (Paik) 

Noise:  0.7  E  Hz‘^^  <  1  Hz  sampling 

- >  All  signals  of  exploration  interest  are  detectable 

A  one  axis  gradiometer  is  sufficient  if  the  signal 

SOURCE  IS  DIRECTLY  BELOW  THE  FLIGHT  PASS.  FULL 
TENSOR  GRADIOMETER  WOULD  INCREASE  CHANCES  FOR 
DETECTION  OF  OFF-PROFILE  SOURCES. 

Future  superconducting 

Noise:  0.07  Hz"^^  <  1  Hz  sampling  rate 

- >  Offers  resolution  of  much  shorter  wavelengths  and 

possibly  new  applications, 

higher  flying  altitude  and  wider  profile  spacing 
possible  for  exploratory  survey. 


CONCLUSIONS 


Speed^  resolution^  and  area  coverage  are  decisive 
ADVANTAGES  OF  AIRBORNE  GRADIOMETRY  OVER  CONVENTIONAL 
GRAVITY  EXPLORATION  METHODS. 

Attenuation  of  gradients  with  height  is  a  problem  for 
SHORT  WAVELENGTH  FEATURES.  A  FLYING  ALTITUDE  OF  200  M 
IS  REQUIRED  TO  RESOLVE  WAVELENGTHS  DOWN  TO  300  M. 

Current  and  near  future  systems  can  resolve  most  signals 

IMPORTANT  IN  PETROLEUM  PROSPECTING  AT  FLYING  ALTITUDES 

OF  200  -  300  M. 

Full  tensor  systems  are  better  than  one  axis  systems 

FOR  INTERPRETING  SPECIFIC  STRUCTURES  AND  FOR  DETECTING 
OFF-PROFILE  SIGNALS. 

Elimination  of  topographical  and  geological  noise  is 
NECESSARY  TO  <  0.5  E  PER  KM  AT  FLYING  ALTITUDE. 


Signal  detection  at  flight  level  is  preferable  to  downward 

CONTINUATION. 

Superconducting  gradiometers  are  advantageous  in  these 

APPLICATIONS  because  OF  HIGH  SAMPLING  RATE  AND  LOWER 


NOISE  LEVEL. 


PAPER  TITLE: 


REQUIREMENTS  FOR  THE  USE  OF  AIRBORNE  GRAVITY  GRADIOMETRY 
IN  GEOPHYSICAL  EXPLORATION 


^^’EAKERS  NAME:  Klaus-Peter  Schwarz 
Questions  and  Comments: 

Anthony  R.  Barringer:  In  airborne  geophysical  surveys  it  is  quite  normal 
to  fly  at  survey  heights  as  low  as  600  ft.  However,  turbulence  increases  a 
great  deal  at  these  heights.  What  are  your  comments? 

SPEAKERS  RESPONSE:  It  may  be  possible  to  stabilize  aircraft  using  wingtip 
devices.  Ihe  problem  is  appreciated  and  there  is  a  trade-off  between  flying 
height  and  turbulence  problems . 

Ho  Jung  Palk:  If  you  have  a  gradiometer  with  higher  sensitivity, 

could  you  take  advantage  of  the  sensitivity  and  fly  at  a  higher  altitude  to 

suppress  the  air  turbulence  noise? 

SPEAKERS  RESPONSE:  Yes,  one  can  always  take  advantage  of  higher  sensitivity. 

Dave  Sonnabend:  What  is  the  effect  of  near  surface  features  causing 

gravity  noise  masking  signals  from  deeper  features  -  maybe  we  want  to  fly  a 
little  higher. 

SPEAKERS  RESPONSE:  Agree  that  one  must  optimize  altitude;  but  even  with  full 
tensor  measured  everywhere,  cannot  sort  this  out  completely. 


QUICK  REVIEW  OF  GRADIOMETER-AIDED  LAND  NAVIGATION 


W.  G.  Heller 

The  Analytic  Sciences  Corporation 
One  Jacob  Way 
Reading,  MA  01867 


ABSTRACT 

Real-time  Inertial  navigation  compensation  for  gravity  disturbances  is  one 
of  the  most  demanding  applications  of  gradiometry.  This  discussion  reviews  the 
potential  contribution  of  existing  gradlometer  Instruments  to  land  navigation, 
namely  the  benefits  which. accrue  to  gravity  vector  estimation  and  inertial  platform 
tilt  determination.  Real-time  gravity  vector  recovery  accuracy  is  seen  to  be 
strongly  dependent  on  the  quality  of  initialization  gravity  data,  inertial  system 
and  gradlometer  accuracy  as  well  as  elapsed  time  from  initialization.  Position 
and  velocity  accuracy  is  observed  not  to  be  significantly  enhanced  with  a 
gradlometer.  The  sharp  contrast  between  real-time  determination  of  the 
gravity  field  and  the  accuracy  obtainable  through  post-traverse  smoothing 
motivates  combining  the  gravity  gradlometer  measurements  with  gravity  vector 
map.  The  sensitivity  of  real-time  gravity  determination  to  gradlometer 
errors  also  motivates  further  development  of  terrestrial  gradlometer  Instruments 
to  achieve  noise  densities  better  than  1.0  E^/Hz. 


THE  ANALYTIC  SCIENCES  CCRPORATION 


SP-4423-30 


QUICK  REVIEW  OF 
GRADIOMETER-AIDED 
LAND  NAVIGATION 


11-12  February  1986 


Prepared  for: 

FOURTEENTH  MOVING  BASE  GRAVITY  GRADIOMETER  REVIEW 
United  States  Air  Force  Academy 

Colorado 


Prepared  by: 
W.G.  Heller 


THE  ANALYTIC  SCIENCES  CORPORATION 
One  Jacob  Way 

Reading,  Massachusetts  01867 


LAND  NAVIGATION  BACKGROUND 
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TYPICAL  REAL-TIME  LAND  INERTIAL  NAVIGATION 
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APPLICATIONS  OF  SUPERCONDUCTING  GRAVITY  GRADIOHETER  SYSTEM 


TOWARD  INERTIAL  GUIDANCE  AND  FUNDAMENTAL  SCIENCE 

H.  A.  Chan 
M .  V .  Moody 
H.  J.  Paik 

Dept  of  Physics  and  Astronomy 
University  of  Maryland 
College  Park,  Maryland  20742 


ABSTRACT 

In  addition  to  the  gravity  survey  application  which  is  the  primary  objective 
^f  the  NASA  program,  the  superconducting  gravity  gradiometer  system  can  be  used 
as  a  gradlometer-aided  inertial  navigation  system  and  as  a  detector  for  a  series 
of  fundamental  physics  experiments  both  in  the  laboratory  and  in  the  earth  orbit. 
These  applications  will  be  discussed. 
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APPLICATIONS  OF  SUPERCONDUCTING  GRAVITY  GRADIOMETF.R  SYSTFM 
TOWARD  INERTIAL  GUIDANCE  AND  FUNDAMENTAL  SCIENCE 


PEAKERS  NAME: 


Ho  Jung  Palk 


Questions  and  Comments: 

r 

Chris  Harrison:  The  Stacey  and  Eotvos  experiments  (as  revised  hy  Fischbach) 

disagree  by  a  factor  of  about  20.  This  is  an  order  of  magnitude  disagreement . 

»  SPEAKERS  RESPONSE:  Speaker  thought  I  (Chris  Harrison)  was  talking  about  the 
Dlcke  experiment. 

Milton  Trageser:  You  showed  a  slide  showing  the  response  to  the  swinging 

pendulum.  Over  how  many  pendulum  swings  was  this  average? 

SPEAKERS  RESPONSE:  500  swings  -  2000  second. 

Paul  Keyser:  With  regard  to  the  inconsistency  between  Dicke  and  Fischback's 

version  of  Eotvos,  it  is  Important  to  realize  that  Fischback's  postulates  a 
short-range  force  so  that  only  ElotvcJs'  results  ar  applicable.  It  is  also 
important  to  note  that  Fischback  left  out  a  numbe .  of  data  points  which  make 
his  curve  less  impressive,  as  I  will  mention  in  my  talk. 

Daniel  DeBra :  How  can  there  be  violation  for  characteristic  lengths  on 

either  side  of  Panov's  data  which  claims  no  violation? 


SPEAKERS  RESPONSE:  Panov's  data  doesn't  overlay  the  length-sensitivity  where 
violation  is  claimed. 


EFFICIENT  GRAVITY  GRADIENT  DATA  GATHERING 


M.  Bllello 
J.  B.  Breakwell 
D.  B.  DeBra 

Department  of  Aeronautics  and  Astronautics 
Stanford  University 
Stanford,  California  94305 


ABSTRACT 

We  are  Interested  In  how  one  can  separate  the  variations  In  gravity  field 
from  the  measurement  noise  when  making  a  survey.  Given  a  survey  pattern  In 
which  the  path  of  the  instrument  crosses  itself  (as  it  does  in  a  series  of 
orthogonal  tracks),  there  are  a  discrete  number  of  instants  at  which  the  meas¬ 
urements  should  be  Identical.  We  have  examined  a  number  of  different  sequences 
in  generating  the  survey  pattern  to  vary  the  times  at  which  these  Identical 
conditions  occur.  The  conjecture  was  that  an  appropriate  choice  of  pattern 
could  take  advantage  of  the  time  characteristics  of  the  measurement  noise  In 
permitting  a  separation  of  noise  from  gravity  data.  We  show  the  results  as  a 
function  of  the  correlation  time  of  the  measurement  noise  for  a  simple  model  of 
the  gravity  field.  For  noise  varying  from  uncorrelated  to  a  correlation  time 
comparable  to  the  survey  time,  the  variation  Is  approximately  10%.  Large 
differences  in  accuracy  of  reconstruction  do  not  appear  likely  since  our  results 
give  variation  between  paths  of  approximately  2%  for  two  very  dissimilar  paths 
through  the  same  grid.  Thus  the  conjecture  has  not  been  borne  out. 


Efficient  Gravity  Gradient  Data  Gathering 


Introdnction 

The  modern  interest  in  measuring  gravity  gradients  began  in  the  late  1950b  motivated 
by  determining  the  vertical  in  a  satellite.  Evly  papers  considering  the  analytical  aspect 
of  gradient  determination  were  followed  in  the  next  decade  by  a  number  of  innovative 
approaches  in  how  such  an  instrument  might  be  built.  The  revolution  in  gradiometry  was 
to  make  the  measurements  in  a  moving  vehicle  and/or  in  a  satellite  without  the  gravity 
needed  for  the  geophysical  pendulum  instruments.  An  instrument  developed  at  the  Bell 
Aerosystems  was  chosen  for  field  application  for  improvements  in  navigation.  This  instru* 
ment  has  been  very  successful  in  its  early  fidd  tests  and  is  in  production  for  deployment. 
As  a  result  of  this  success  for  the  navigation  mission,  the  Defense  Mapping  Agency  (DMA) 
through  the  Air  Force  Geophysical  laboratory  (AFGL)  began  the  modification  of  this  in* 
strument  for  gravity  gradient  measuremoits  for  gravity  survey  work.  Many  people  have 
subsequently  contributed  to  the  development  of  a  survey  plan  and  techniques  for  utilisa¬ 
tion  of  such  an  instrument.  In  this  paper  we  explore  the  possibility  that  pven  an  area 
to  be  surveyed  and  a  track  spacing  that  has  been  determine  by  the  necessary  resolution 
of  gravity  data,  there  might  be  improvements  in  accuracy  depending  upon  the  form  of  a 
grid  pattern  used  in  overflying  the  area.  The  conjecture  is  based  on  the  fact  that  instru¬ 
ment  noise,  whether  described  in  the  time  domain  or  spectrally,  may  be  different  than  the 
equivalent  noise  associated  with  gravity  fields  for  a  pven  velocity  of  the  vehicle  during  the 
survey.  When  a  survey  is  performed  with  a  grid  in  which  tracks  cross  each  other,  there 
are  a  discrete  number  of  crossings  at  which  the  measurements  should  be  the  same  in  both 
directions.  Different  patterns  provide  a  different  distribution  in  time  of  when  these  points 
of  identical  measurement  occur.  It  is  this  variation  In  the  distribution  and  time  which 
could  make  a  difference  in  b«ng  able  to  separate  signal  from  noise. 


Models 

As  indicated  in  the  introduction,  the  spectral  characteristics  of  the  gravity  field  and 
of  the  instrument  will  have  an  influence  on  the  separability  of  the  gravity  information 
from  the  instrument  noise.  With  the  amount  of  experimental  data  that  exists  from  the 
laboratory  and  early  field  trials,  it  would  be  possible  to  give  a  good  model  of  the  expected 
noise  from  a  gravity  gradiometer.  However,  to  Investigate  the  potential  for  improvement 
one  can  start  with  a  much  simpler  model  of  the  instrument  noise  and  vary  its  parameters 
to  see  whether  or  not  significant  improvements  are  possible.  We  have  chosen  the  latter  ap¬ 
proach  to  investigate  the  feasibility  of  improvement  with  the  expectation  that  if  significant 
improvements  appear  possible  we  would  then  improve  the  model  using  available  empirical 
data. 


Spectral  Cliaracteristics  of  the  Field 

We  have  used  a  model  of  the  gravity  gradient  field  that  allows  us  to  determine  the 
spatial  correlations  of  the  gravity  gradient.  (J.V.  Breakwell  (l)). 
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Using  an  appraximation  of  flat  earth,  we  can  write: 

U{(3th)  —  e~^U(£3to)  where  U((3,o)  is  the  Fonrier  transform  of  (/(x,y,o),  potential 
on  the  reference  surface  of  the  earth  and  U(J,h)  is  the  Fourier  transform  of  U{x,y,h), 
gravity  potential  at  adtitude  h. 

Then  the  gravity  gradients  components  are  pven 


where 


e 


-hu 


-ju,u 


U{u,o) 


U  =  (u„Uf) 

U  =  yjul  +  W* 


(1) 


FVom  Heller’s  model  referenced  in  [1],  we  get  the  spectral  density  of  U{x,y,o)  with 
correlation  distance  D,  : 


^u.(w)  =  ^u,W)  — 


«=i 


Equation  (1)  can  be  viewed  as  a  representation  of  a  linear  system  with  U{u},o)  as 
input  and 


-;w.w 

L-M"J 


as  the  transfer  function. 

Then  we  can  compute  the  spectral  densities  of  the  gravity  gradient  components  at 
altitude  h: 


^Uo{u) 


3 


By  taking  the  inverse  Fourier  transform,  we  can  determine  the  auto*  correlation  func* 
tions  for  the  gradients,  say  Sv{x,  y,  h). 

Example:  Say  we  want  to  compute  Su^^{x,y,h) 
we  have 


s 

»=1 


then 


St...(*.».*)  =  E  *ij  du.iu, 

,'st  “* 


or 


Su..(r,6,h)  =  f^  duda 

,_J  /o  70 


that  is 


.*1 

In  the  special  case  of  a  flight  path  over  a  point  grid,  we  need  to  compute  Sp  s  £[s,sj], 

r*Fi 

where  S,  is  the  sequence  of  signals  we  want  to  estimate  Sp  = 


Let’s  suppose  that  we  are  measuring  the  component  U,,  of  the  gradient  then: 


where 

ro-  =  ll^ll, 

is  the  distance  between  points  Pi  and  P^. 


Gravity  Survey 


To  perform  a  gravity  survey,  the  craft  which  carries  the  instruments  follows  a  partic¬ 
ular  path.  In  the  simple  case  of  a  square  survey  area,  a  possible  strategy  is  to  fly  parallel 
tracks  as  shown  in  Figure  1. 


However,  in  order  to  remove  drifts  and  red  noise  from  the  measurements,  a  better 
way  is  to  make  cross  checks,  taking  two  measurements  at  two  different  times  at  the  same 
point.  The  grid  of  Figure  2.1  is  an  example  of  this  type  of  flight.  Also  shown  is  the  time 
of  second  crossing.  Figure  2.2,  and  the  time  between  the  two  crossings  versus  the  point 
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of  interest,  Figure  2.3.  One  can  see  that  for  the  path  of  Figure  2.2,  the  crossings  occur 
essentially  during  the  second  half  of  the  total  survey  time  T  and  that  when  they  begin  to 
occur,  it  is  in  such  a  way  that  they  are  close  to  each  other  in  space. 

In  order  to  get  a  better  time  and  space  distribution  of  second  crossings,  a  path  such 
as  the  one  shown  in  Figure  3.1  might  be  of  interest.  Here,  a  row  or  column  is  dipped  at 
each  pass,  and  the  effect  can  be  seen  in  Figures  3.2  and  3.3.  Basically,  second  crossings 
occur  earlier  and  two  consecutive  ones  are  more  likely  to  be  spread  in  time.  Another 
advantage  of  this  kind  of  path  is  the  possibility  to  continue  to  make  measurements  while 
turning  between  two  tracks.  If  one  row  or  column  (or  more)  is  skipped,  then  the  radius  of 
curvature  in  the  turning  is  bigger,  so  that  both  the  bank  angle  and  the  induced  acceleration 
are  smaller.  This  may  allow  the  instrument  platform  to  remain  in  tolerable  perturbations 
and  compensations  may  be  possible. 

In  view  of  the  disappointing  results  that  we  are  about  to  give,  we  did  not  pursue  the 
question  of  efficiency  due  to  variations  in  the  radius  of  turns,  nor  did  we  carry  the  study 
to  include  the  effect  of  mass  attraction  and  error  modeling  on  the  instrument. 

Criteria  for  Comparison 

Our  purpose  is  to  get  an  estimate  of  the  gravity  gradient  at  the  grid  points  with  the 
smallest  error-standard  deviation.  Since  all  points  are  a  priori  of  equal  importance,  we 
take  as  the  performance  criterion  the  arithmetic  mean  of  the  standard  deviation  obtained 
at  each  point,  that  is: 


'ptt 


1  ^ 
"  i=i 


where  (r,-  =  ^/Pi  and  Pi  is  the  variance  of  the  error  in  the  gravity  quantity  at  point 
i:  Pi  =  E[{$f,  -  N  is  the  number  of  points  on  the  grid.  Thus,  we  will  be  considering 
as  the  best  path  the  one  that  minimises  the  criterion 


Theory 

The  gradiometer  output  signals  consist  of  the  sum  of  a  signal  to  be  estimated  (gravity 
quantity)  and  the  noise  inherent  in  the  instrument. 

y  =  3  +  n  where  s  is  any  one  of  the  gravity  gradient  components  and  n  is  the 
instrument  noise.  If  we  take  M  measurements  at  A#  different  times,  we  have  in  vector 
form: 


y  =  s  +  n  where  y  = 

yi 

s  = 

... 

;n  = 

ni 

.Vm. 

.»»w. 

where  n,  is  the  instrument  noise  at  time  etc.  ... 

If  the  pattern  is  a  square  grid  with  intersecting  points,  then  M  =  2p*  where  p  is  the 
number  of  points  on  the  side  of  the  square  grid. 


5 


The  aketehee  above  show  4  bjr  4  grids.  The  speed  of  the 
craft  la  uniform  and  the  turning  times  are  negleeted.  la  Fig. 
2.2  and  3.2,  the  time  of  the  second  crossing  at  each  point 
(from  1  to  16)  is  plotted,  while  Pigs.  2.S  and  24  show  the 
time  between  the  two  crossings  for  each  point  (ftrom  1  to 


We  assume  for  simplicity  a  linear  estimate  from  the  observations: 

3  =  Ky  where  K  u  9Sk  M  x  M  gain  matrix.  This  is  a  smoothing  formula  where  we 
use  all  the  collected  data  to  estimate  the  gravity  quantity  at  each  point. 

The  error  in  the  estimate  isi  =  s-3or3  =  {/  —  K)3  -  Kn. 

Then  the  covariance  matrix  of  the  error,  say  P,  can  be  computed: 

P  =  £?I33’’l  =  (/  -  K)S{I  -  Kf  +  KNK^  -(I-  K)Elsn^]K^  -  KElns^](I  -  K^) 
where  5  i  £[««»■)  and  JV  4  £?(nn’’|. 

The  gravity  signal  5  and  the  instrument  noise  being  uncorrelated,  the  formula  reduces 

to: 


P  =  (/  -  K)SiI  -  Kf  + 


Then  we  choose  the  gain  matrix  K  that  minimises  the  trace  of  P  (least  squares 
estimate)  that  is: 


d{tr{P))  =  tr((-(/  -  K)S  +  KN)dK'^  +  dKi-S{I  -  Kf  +  =  0 

This  yields  -(/  -•  K)S  +  KN  =  0  or  if  =  S{S  +  N)"!  whenever  the  inverse  exists 
and  in  this  case  K  exists  and  is  unique.  We  remark  that  if  (5  +  iV)  is  non  invertible  then 
Pjyy^l  is  non  invertible.  Minimising  every  term  Pa  leads  to  the  same  gain  matrix  K.  The 
linear  least  squares  estimate  is  then  deduced  3  =  5(5  +  N)~'y.  The  performance  of  the 
estimate  is  judged  upon  the  error  covariance  matrix  and  more  precisely  on  the  diagonal 
entries  of  this  matrix.  Substituting  for  K  in  the  expression  of  P,  we  get: 

P  =  5(5  +  iV)-*jV 

In  addition  to  the  fact  that  n  and  s  are  uncorrelated,  we  have  implicitly  assumed  that 
s  and  n  are  sero-mean  signals.  If  this  is  not  the  case,  (£(s)  0  and/or  P(n]  0  but  still 

P[ns^]  =  P[sn^]  =  0),  then  the  formulae  are  modified  in  such  a  way  that  we  replace  the 
random  variables  with  their  centered  counterparts,  namely; 

where  (y  =  s  +  n) 


(  3  =  £(s)  +  if(y  - -B(y)) 

I  K  =  5‘(5*  +  NT' 

{  P  =  5*(5*  +  NT'ff' 


with 

5*4£[m’’)  -  E\8]E\a‘^] 


e 


JVi£(nn’'l  -  E[n\E\n^] 


Then  for  a  particular  pattern  that  links  times  to  points,  we  associate  the  variance 
with  the  point  which  is  flown  over  at  time  f/. 

However,  for  a  grid  with  crossed  points,  it  turns  ont  that  it  is  never  necessary  to  take 
the  inverse  of  the  M  x  M  matrix  (5  +  N)  because  as  can  be  expected,  there  are  a  lot  of 
redundancies  in  the  matrix  P  computed  as  P  —  S(S  +  For  example,  if  at  times  j 

and  k  the  same  point  is  flown  over  (with  j  ^  k),  then  obviously  P(t^-,t|)  =s  P(<ft,ti)  V/;  in 
particular,  P(ti,ty)  =  P(t*,t*). 

We  detail  this  in  the  next  section  on  preliminary  numerical  results. 

Numerical  Results 

We  take  for  our  example  p  =  4  and  there  are  16  =  points  on  the  grid  and  we  show 
first  how  to  reduce  the  sise  of  the  matrix  to  be  inverted  (5  +  N)  (the  path  lasts  M  units 
of  time). 

Let 


■ 

and  S,  = 

• 

• 

• 

also  n,  s 

• 

» 

‘n,i' 

• 

• 

.**«*,. 

where  the  subscripts  t  stand  for  time  and  p  for  points  (N  s  p^) 
then 


{ 


si  =  Fa, 
n,  =  Fn, 


where  P  is  the  it/  x  p’  matrix  that  maps  the  points  to  the  times,  i.e.,  F(i,j)  =  1  if 
point  j  is  flown  over  at  time  and  0  otherwise. 

F  is  full  rank  and  let  Fi  be  the  pseudo-inverse  of  F(Fi  is  p^  x  Af)  then  we  can  write 

/  5,  =  FS,F^  I  S,  =  P,5,f7 

( iV,  =  Ffr,f^  \  N,  =  Fin,f7 


where 

5  =  P(ss*’l 
N  s=  Plnn**) 


FVom  previous  results  we  had: 


P,  =  5,(5,  +  N,)-^Nt 
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which  yield! 

P,  =  FS,F^[FS,f^  +  FN,F^]-^FN,F'^ 
or 

P,  =  F5,F^1P(S,  +  N,)F'^]-^FN,F'^ 
but 

f^[F{S,  +  N,)F^\-'F  =  (5,  + a;)-'. 

Then 


Pi  =  FS,{S,  +  N,)''N,F^  =P  P^P^  where  P,  ^  5,(5,  +  Ap)-‘Ap. 

is  a  X  matrix  the  diagonal  entries  of  which  are  repeated  in  the  diagonal  of  Pi. 
P,  gives  directly  the  covariance  of  the  gravity  gradient  at  the  points  of  interest. 

For  the  numerical  example,  we  chose  a  4  x  4  grid  with  two  different  paths  and  we 
wish  to  compare  the  performances  using  the  criterion  mentioned  earlier.  We  have  first  to 
define  the  covariance  matrices  N  and  5  and  to  construct  the  F  matrix  for  the  two  different 
paths. 

The  models  used  for  the  random  signals  n  and  s  are  exponentially  correlated.  That  is, 
the  entries  of  the  covariance  matrix  Ni  vary  as  the  exponential  of  the  time  difference  and 
the  entries  of  the  covariance  matrix  5,  vary  as  the  exponential  of  the  distance,  namely: 

ir,/| 

=  e  '  and  5,(»,y)  =  e~^ 

where  r  and  are  correlation  time  and  correlation  distance,  respectively. 

These  models  do  not  claim  to  be  accurate  but  represent  only  a  first  try  to  get  numerical 
performance. 

Then  we  compute  P,  =  5,(5,  +  A,)“*A,  to  determine  the  variance  of  the  error 
associated  with  the  gravity  gradient  at  each  point  of  the  grid. 

For  the  two  paths,  we  plot  the  mean  of  the  standard  deviation  versus  r(r  =  0  corre¬ 
sponds  to  a  white  noise). 

ConcliuioB: 

The  spectral  models  of  instrument  noise  and  gravity  gradient  signal  we  used  in  our 
simulations  may  not  be  realistic  and  this  marks  the  limitation  of  the  results  we  got.  How¬ 
ever,  in  the  special  case  of  exponential  correlated  signals,  they  allow  us  to  answer  the 
question  of  the  best  path  (among  specified  ones)  according  to  the  criterion  we  defined.  In 
terms  of  times  of  second  crossing  and  times  between  crossings,  the  two  paths  chosen  for 
the  simulation  can  be  described  as  Very*  different.  Surprisingly  enough,  the  performances 
for  the  two  paths  are  close  to  each  other  for  the  range  of  correlation  times  we  have  run. 


s 


However,  the  gap  it  getting  wider  in  favor  of  path  1  when  the  correlation  time  gets  larger 
but  the  performance  of  path  1  is  only  1.5  %  better  for  r  13  unite  of  time.*  In  these 
conditions,  the  choice  of  a  "better*  path  appears  not  to  be  an  issue. 

Our  final  remaric  concerns  the  nature  of  the  instrument  noise.  The  way  it  has  been 
modelled  assumed  that  it  was  stationary  (in  particular  constant  variance  at  any  time);  if 
this  is  not  the  case,  quite  different  results  may  occur;  for  example,  the  importance  of  early 
crossings  increases. 


*  1  unit  of  time  is  the  time  required  to  fly  from  a  pout  to  the  next  one. 


» 


lieaN  RMS  versus  Cor  re  1  a-fc J on  Time  of  iNs-fc 


f'Joi  Ml  y  Jsn-aa  ^a^l^ 


J. 


Ill  Breakwell,  J.V.,  *Satellite  Delerminatioa  of  Short  Wavel«tgth  Grarity  Variatioiif 
of  the  Astronaiaual  Seience$,  No.  4,  pp.  32d-344,  Oct-Dec.  1®79.  ’ 
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PAPER  TITLE: 


EFFICIENT  GRAVITY  GRADIENT  DATA  GATHERING 


SPEAKERS  NAME:  Michel  Bllello 
[Questions  and  Comments: 

Peter  Uglnclus:  !•  If  the  noise  were  not  correlated,  would  there  be  any 
difference  for  different  survey  patterns?  2.  What  determines  the  position 
of  the  maximum  in  the  plot  of  error  covariance  vs.  noise  correlation  time? 

SPEAKERS  RESPONSE;  1.  No.  2.  The  ratio  of  correlation  times  for  noise 
and  signal. 
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USE  OF  TERRAIN  ELEVATION  DATA  IN  AIRBORNE  GRADIOMETRY 

J.  L.  Center,  Jr. 

Geospace  Systems  Corporation 
Wellesley  Office  Park 
40  William  Street 
Wellesley,  Mass.  02181 


ABSTRACT 


Terrai 
This  paper 
Suggests  a 


n  elevation  data  is  needed  to  relate  airborne  survey  data  to  ground 
quantifies  the  errors  introduced  if  terrain  effects  are  ignored.  It 
method  for  operational  use  of  terrain  elevation  data. 


data. 

also 


coM^OMn 


WHY  ARE  TERRAIN  EFFECTS  IMPORTANT? 
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SURVEY  ERROR  ANALYSIS  FLOW 


PAPER  TITLE: 


USE  OF  TERRAIN  ELEVATION  DATA  IN  AIRBORNE  GRADIOMETRY 


SPEAKERS  NAME:  Julian  Center 


estions  and  Comments: 


^  Dan  DeBra  :  Is  a  single  site  an  adequate  basis  for  verifying  a  GGSS  and 
data  reduction  systems  given  the  disagreement  based  on  simulation? 


Terry  Fundak:  The  site  chosen,  Bakers'  Peak,  has  a  variety  of  terrain- 
some  with  lots  of  topographic  relief  and  relatively  little  gravity  structure 
and  some  with  flat  terrain  which  Is  geologically  active.  White  Sands  or 
another  site  is  possibility. 


Chris  Jekell:  How  do  you  know  how  much  terrain  correction  to  subtract  from  the 
data  without  compromising  the  Instrument  test  In  terms  of  what  It  will  sense. 

SPEAKERS  RESPONSE:  There  should  be  no  problem.  If  you  subtract  out  a  major 
portion  of  the  signal  and  you  have  no  problems  In  data  reduction  you  know  the 
Instrument  works  (has  sensed  the  field  you  subtracted). 
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THE  GRAVIIY  GRADIOMETER  SURVEY  SYSTEM 

T.  J.  Fundak 

Air  Force  Geophysics  Laboratory 
Earth  Sciences  Division 
Hanscom  AFB,  MA  01731 


ABSTRACT 

The  Gravity  Gradiometer  Survey  System  (GGSS)  is  now  entering  its  first 
phase  of  testing  as  a  system.  Over  the  next  year  a  series  of  tests  will  be 
conducted  to  understand  the  limits  and  characteristics  of  the  instrument  and 
^he  gradient  environment  in  which  it  operates.  This  presentation  will  outline 
the  status  of  the  GGSS  development  program  and  highlight  the  events  to  take 
place  over  the  next  year  in  the  GGSS  Test  Program. 
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The  Gravity  Gradiometer  Survey  System  (GGSS)  is  a  self-contained  gravity 
measuring  system  designed  to  provide  gravity  gradients  from  a  moving  vehicle. 
The  GGSS  combines  state  of  the  art  inertial  technology  and  satellite  based 
navigation  to  produce  a  system  vrlth  the  capability  to  provide  gravity  data  in 
all  types  of  environments  throughout  the  world.  The  GGSS  is  being  developed  to 
perform  gravity  surveys  from  an  airborne  platform  or  in  a  land  vehicle  with  a 
recovery  of  the  gravity  vector  to  an  accuracy  of  better  than  1.0  milligals 
magnitude  and  .18  arcsecs  in  each  deflection  component.  For  more  details  about 
the  GGSS  specifications  see  Table  1. 


A  gravity  gradiometer  works  by  differencing  the  applied  forces  (both  inertial 
and  gravitational)  across  the  distance  between  two  acceleration  measurements. 

Most  of  the  linear  inertial  acceleration  environment  is  naturally  rejected  by 
this  method,  leaving  only  the  difference  in  gravitational  force  across  that 
distance:  the  gravity  gradient  (Figure  1).  By  properly  Integrating  the  gravity 
gradient  along  a  survey  track,  the  relative  value  of  gravity  between  two  points 
can  be  found.  In  contrast  to  a  gravimeter,  however,  all  three  components  of 
the  gravity  vector  can  be  found  with  a  tensor  gradiometer,  such  as  the  GGSS. 


Like  all  geophysical  quantities,  a  measurement  is  only  valuable  when  the 
location  of  the  measurement  is  known.  In  the  case  of  a  gravimeter,  the  estimate 
of  gravity  at  the  1  milllgal  accuracy  level  is  only  valid  and  useful  if  the 
location  of  the  measurement  is  known  to  about  a  meter  in  vertical  and  about  a 
hundred  meters  in  horizontal  position.  However,  in  contrast  to  gravimeter 
measurements,  in  order  to  produce  similar  accuracy  in  the  recovery  of  gravity 
the  gravity  gradients  are  needed  to  only  one  part  in  3000.  In  an  airborne 
environment,  this  means  that  the  need  for  positioning  a  measurement  of  gravity 
gradients  in  space  is  not  nearly  as  stringent  as  for  a  conventional  gravimeter 
survey.  One  hundred  meters  of  absolute  position  error  is  acceptable  to  recover 
gravity  to  better  than  one  milllgal  in  accuracy.  A  satellite  based  navigation 
system,  such  as  Global  Positioning  System  (GPS),  can  routinely  provide  position 
information  with  this  level  of  accuracy. 


Near  the  surface  of  the  earth,  however,  positioning  becomes  more  critical. 
While  the  gravity  field  is  relatively  well  characterized  even  at  the  high 
frequencies,  the  high  frequency  variations  of  the  gradient  field  at  the  1 
Hot VOS  level  are  not  well  characterized  to  date.  Figure  2  indicates  that  many 
common  objects  produce  significant  signals  that  can  normally  be  Ignored  in 
measuring  the  gravity  field.  It  is  expected,  however,  that  the  power  spectrum 
of  the  horizontal  components  of  the  gradient  field  will  not  have  significant 
power  at  wavelengths  of  tens  of  meters.  To  provide  the  meter  level  navigation 
accuracy  that  may  be  needed  will  require  a  supplement  to  the  GPS  system. 


HAROWARE 


The  heart  of  the  GGSS  system  consists  of  a  gyro-stabilized  platform  with 
the  three  orthogonal  gravity  gradient  Instruments  mounted  in  an  umbrella 
configuration.  Six  electronics  racks  containing  an  airborne  computer,  a  GPS 
Navigation  System,  power  supplies,  control  electronics,  a  vehicle  interface 
buffer,  and  data  recording  equipment  service  the  gradiometer  and  platform, 
and  record  data  for  post  mission  reduction. 


* 


1  Eotvos 


ten  to  the  minus  nine  per  second  sqiiared 
0.1  milllgal  over  1  kilometer. 
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The  Gravity  Gradioaeter  Instrument  and  the  Platform 

The  basic  gravity  gradiometer  instruaent  consists  of  2  pairs  of  specially 
selected  accelerometers  (figure  3).  A  pair  of  accelerometers  provides  a 
measurement  of  the  average  gravity  gradient  field  across  the  distance  that 
separates  them.  In  a  dynamic  envlroment,  it  is  necessary  to  use  a  slow  rotation 
to  modulate  the  gravity  gradient  signal.  This  technique  allows  the  use  of 
frequency  domain  signal  extraction  techniques,  so  that  the  signal  to  noise 
ratio  can  be  Improved.  In  this  way,  gradient  signals  many  millions  of  times 
smaller  than  the  acceleration  which  produce  them  can  be  extracted.  This  method 
has  proven  very  effective  in  the  Navy's  Gravity  Sensing  System  (GSS)  program 
where  the  Bell  Aerospace  gradiometer  is  being  used  to  correct  ship  inertial 
navigation  systems  for  gravity  anomaly  induced  errors. 

Two  sets  of  accelerometer  pairs  are  mounted  as  orthogonal  pairs  at  the  edge 
of  a  rotating  disk  (figure  4).  Theses,  plus  control  electronics,  are  then  mounted 
in  two  half  spherical  shells.  This  makes  up  one  Gravity  Gradiometer  Instrument 
(GGI)  (figure  5)  .  A  set  of  three  GGIs  is  mounted  on  the  gyro-stablized 
platform  mechanized  to  provide  a  north,  east,  down  (NED)  or  "Local  Level" 
environment.  The  platform  maintains  the  local  level  attitude  with  the  aid  of 
position  information  from  the  GPS  receiver  for  the  airborne  survey  and  a  "fifth 
wheel"/lnertlal  mechanization  for  the  land  survey  .  The  stabilized  platform 
itself  acts  as  a  short  term  navigation  system  in  a  free  inertial  mode,  if  the 
GPS  signal  is  lost  for  a  short  period  of  time.  With  this  mechanization,  the 
long  term  navigation  error  (drift)  is  handled  by  the  GPS  system  and  the  short 
term  high  accuracy  navigation  is  handled  by  the  inertial  platform  (two 
two-degree-of-freedom  gyros  and  a  triad  of  accelerometers). 

Electronics  &  Interfaces 

While  the  heart  of  the  GGSS  is  mechanical,  the  brains  are  in  the  electronics. 
Figure  6  shows  a  block  diagram  of  the  major  components  of  the  system  and  the 
signal  flow.  At  the  center  of  all  the  electronics  is  the  ROLM  MSE-14  airborne 
computer.  This  computer  is  used  for  everything  from  platform  control  and 
aircraft  survey  mission  control  to  data  recording  and  system  health  and  control 
functions  monitoring. 

A  majority  of  the  GGSS  electronics  is  housed  in  five  electronics  racks. 

Rack  1  takes  care  of  the  GGI's  needs.  Rack  2  monitors  and  controls  platform 
parameters,  takes  care  of  housekeeping  and  serves  as  an  interface  between  the 
platform,  the  GGIs  and  the  computer.  The  other  three  racks  of  e'.ectronlcs 
contain  the  computer,  a  GPS  reciever,  an  atomic  clock,  data  recording  equipment, 
and  interfaces  to  the  various  aircraft  sensors  and  systems. 

All  of  the  various  interfaces  necessary  to  monitor  the  aircraft  sensors  and 
systems,  and  to  feed  control  signals  to  the  autopilot  and  to  receive  power  for 
the  GGSS  are  funnelled  through  a  vehicle  interface  buffer  (VIB).  The  antenna 
lead  for  the  GPS  reciever,  for  example,  is  interfaced  through  the  VIB  to  the 
Texas  Instruments  GPS  reciever  which  in  turn  is  Interfaced  with  the  computer. 

The  GPS  system  consists  of  a  TI-4I00  GPS  Navigator  and  a  CHU  aircraft  antenna. 

Another  peripheral,  but  significant,  component  is  the  Uninterruptable  Power 
Supply  (UPS)  and  battery  backup  which  provides  power  conditioning  and  emergency 
power  for  the  GGSS.  (The  aircraft  power  environment  is  very  noisy  and  the  GGSS 
is  very  unforgiving  when  it  comes  to  power  surges,  transients,  drop  outs,  etc.) 
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Land  Vehicle  and  Aircraft 


The  GGSS  will  be  Installed  in  a  modified  motorhome  to  perform  the  land  test 
and  in  a  commercial  C-130  Hercules  aircraft  for  the  airborne  survey.  The 
choice  of  aircraft  for  the  GGSS  was  made  to  maximize  testing  and  operational 
flexibility.  The  C-130  aircraft  is  also  well  suited  for  low  level  airborne 
surveys  that  will  be  required  of  the  GGSS  and  will  require  little  modification. 
Aircraft  modification  can  be  a  very  long  and  costly  process  due  to  the  complex 
nature  of  the  aircraft  environment  -  safety  is  always  the  first  priority. 

With  the  GGSS  program,  an  air  transportable  land  vehicle  will  be  loaded  into 
the  C-130  aircraft,  strapped  down,  connected  to  the  various  interface  cabling, 
and  be  flown  on  the  day's  mission.  This  arrangement  should  help  minimize 
"down  time"  of  the  GGSS  due  to  aircraft  related  maintenance  and  repair.  A  back 
up  aircraft  can  easily  be  substituted,  if  extended  maintenance  time  is  required. 
This  is  in  contrast  to  a  system  that  is  “hardwired"  to  an  aircraft. 

SOFTWARE 


The  GGSS  contains  three  levels  of  software:  1)  real-time,  and  two  stages  of 
post  mission  software,  2)  Stage  1  -  time  domain  processing,  and  3)  Stage  2-  space 
domain  processing.  A  functional  block  diagram  of  the  software  is  shown  in  figure 
7.  All  data  reduction  is  accomplished  off-line  with  the  aid  of  two  stages  of 
reduction  software.  Stage  1  Processing  is  common  to  both  the  land  and  airborne 
tests,  whereas  Stage  2  Processing  is  tailored  to  the  particular  survey  vehicle 
and  conditions. 

Real  Time  Software 


The  real  time  software  is  run  on  imbedded  microprocessors  and  a  ROLM  MSE-14 
airborne  computer.  During  the  survey  the  MSE-14  monitors  system  performance, 
controls  the  vehicle  and  platform,  and  records  the  system  outputs.  GGI  and 
other  Instrument  outputs  are  corrected  for  environmental  sensitivities,  such  as 
pressure  and  temperature,  and  are  compensated  for  instrument  bias  and  drift,  in 
real  time.  The  GGSS  Instrument  outputs  are  then  stored  on  magnetic  tape  for 
post  mission  analysis  and  reduction.  The  outputs  of  the  real  time  software  are 
modulated  gradients,  accelerometer  outputs,  platform  angles  and  rates,  GGI 
wheel  angles,  position  information  (latitude,  longitude,  and  altitude),  time, 
and  other  sampled  Instrument  outputs  from  various  aircraft  and  GGSS  sensors. 

Stage  1  Data  Processing 


Stage  1  software  takes  the  modulated  gradiometer  outputs,  platform  parameters, 
and  position  information  and  transforms  these  time  domain  signals  into  estimates 
of  the  gradient  in  the  instrument  coordinates  at  the  location  of  the  measurement 
in  space.  During  this  stage,  compensations  are  made  for  the  self  gradient  of 
the  platform,  static  and  time  varying  vehicle  Induced  gradients,  angular  acceleration 
sensitivity  and  other  residual  instrument  sensitivities.  The  outputs  of  this 
stage  of  processing  are  estimates  of  the  along  track  gradients  and  their  location 
in  space. 

Stage  2  Data  Processing 


During  Stage  2  Data  Processing,  the  output  of  Stage  1  is  transformed  into 
estimates  of  gravity.  Stage  2  Land  software  takes  the  along  track  gradient 
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estimates  and  integrates  them  from  a  “truth  data  point”  to  the  point  where  a 
comparison  is  desired.  Stage  2  Airborne  software  contains  a  set  of  routines 
which  takes  the  gradient  estimates  at  survey  altitude,  grids  them,  compensates 
them  for  track  crossing  mismatch,  incorporates  the  long  wavelength  Information 
contained  in  the  astro-geodetic  tie  points,  and  dotmward  continues  the  data  to 
the  surface  of  the  earth.  Estimates  of  the  gravity  field  are  then  made  for  the 
area  of  overflight.  The  outputs  are  an  estimates  of  the  gravity  vector  on  the 
Earth's  surface. 


TESTING 


Testing  will  be  conducted  in  two  geographic  areas  of  the  United  States. 

The  first  site  was  chosen  for  its  proximity  to  the  Bell  Aerospace  plant.  It 
will  include  the  areas  adjacent  to  Niagara  Falls,  New  York.  The  second  site 
was  picked  for  its  excellent  truth  data  coverage,  limited  topographic  relief, 
nearby  contrasting  areas  of  active  and  benign  gravity  field,  and  good  weather 
and  aviation  environment.  It  is  located  in  the  Texas /Oklahoma  Panhandle  region 
of  the  Central  United  States  (see  figures  8  thru  10). 

The  GGSS  test  program  is  divided  into  three  parts:  Laboratory  testing, 
and  two  "Phases"  of  Land  and  Airborne  Testing.  As  a  matter  of  course,  each 
test  will  build  on  the  knowledge  gained  in  the  previous  one.  Each  will  provide 
a  different  look  at  the  use  of  a  gradlometer,  its  capabilities  and  shortfalls 
and  each  will  provide  unique  challenges  for  testing.  For  example,  the  airborne 
environment  provides  for  an  attenuation  of  the  high  frequency  gravity  gradient 
field,  relatively  low  special  sampling  rates  because  of  the  aircraft  velocity 
and  a  potentially  hostile  acceleration  environment.  The  signal  in  the  land 
environment,  on  the  other  hand,  may  be  dominated  by  the  high  frequency  signal. 
Even  at  relatively  high  sampling  rates  aliasing  may  occur,  thus  requiring 
faster  sampling  rates,  a  slow  vehicle  survey  speed  and  more  accurate  positioning 
information. 

The  way  in  which  the  data  from  each  phase  of  testing  is  reduced  is  also 
significantly  different.  In  the  case  of  the  land  testing,  the  gradients  are 
integrated  between  a  start  and  an  end  point.  Any  error  between  the  known  relative 
values  and  gradlometer  measured  values,  is  a  direct  indication  of  Instrument 
performance.  In  the  airborne  test,  a  grid  of  data  at  flight  altitude  is  used 
to  make  estimates  of  gravity  values.  Downward  continuation  and  estimation 
procedure  errors  may  dominate  the  errors  sources  here.  Spatial  aliasing  errors 
will  also  limit  the  system's  level  of  performance.  The  aircraft's  relatively 
high  velocity  gives  a  relatively  low  spatial  sampling  rate.  This,  along  with 
the  finite  special  sampling  of  the  grid  inherently  limit  the  spatial  resolution. 

The  ability  to  optimally  reduce  the  gradlometer  data  post-mission,  while  an 
important  question  for  comparison  to  truth  data,  is  not  a  critical  factor  when 
only  GGI  noise  performance  is  being  gauged.  Therefore,  two  types  of  tests  are 
performed  for  each  of  the  two  testing  environments  (land  and  airborne).  The 
first  test  type,  "Phase  I  testing",  gauges  the  system's  sensitivity  to  various 
vehicle  and  environmental  conditions,  determining  the  true  instrument  noise 
performance.  The  second  test  type  "Phase  II  testing”  is  designed  to  assess  the 
performance  of  this  survey  tool  under  a  nominal  set  of  "operational  conditions". 

Laboratory  Testing 

Laboratory  Testing  consists  of  a  series  of  tests  which  are  designed  to  take 
the  individual  components  of  the  GGSS  through  a  series  of  calibration  and 
performance  verification  runs.  Component  level  testing  is  followed  by  the 
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Integration  and  systems  level  testing.  Static  and  Dynamic  Tests  are  performed 
with  the  aid  of  a  Scorsby  Table  and  various  laboratory  calibration  and  checkout 
equipment.  The  various  operating  modes  of  the  system  are  checked  and  a  set  of 
error  coefficients  for  the  platform  Induced  self  gradients  is  produced,  in 
addition  to  the  bias  and  trend  e?itlmates. 


Phase  I  Testing 


Phase  I  is  both  an  Instrument  performance  test  and  an  optimization  tool  for 
Phase  II  testing.  It  will  allow  for  more  effective  "survey  planning".  During 
this  phase  of  testing,  each  of  the  controllable  parameters  is  varied  and  optimized 
for  use  in  Phase  II  testing.  Each  day's  run  is  used  to  estimate  the  system's 
expected  performance  and  sensitivities  so  that  a  simulated  Phase  II  test  performance 
result  can  be  estimated.  An  understanding  of  the  characteristics  of  the  land 
and  airborne  gradient  environments  will  also  be  gained  during  this  phase  of 
testing.  These  characteristics  are  also  fed  into  the  simulation  to  understand 
how  the  gradlometer  will  be  affected  by  differing  gravity  field  power  spectrums. 

The  software  Is  then  optimized  for  Phase  II  testing  and  operations.  In  summary. 
Phase  I  will  evaluate  the  conditions  under  which  the  Instrument  has  ability  to 
give  consistent  results,  l.e.  can  maintain  the  noise  at  the  requisite  level. 

The  testing  strategy  for  this  phase  of  testing  can  be  summed  up  in  two 
words  -  repeat  tracks.  To  evaluate  the  instrument  noise  characteristics,  the 
GGSS  is  required  to  traverse  a  particular  track  many  times.  Statistical  analysis 
will  then  be  applied  to  the  data  to  determine  the  instrument  noise  floor,  drift 
patterns  and  environmental  sensitivities. 


Phase  I  land  testing  will  be  conducted  on  a  set  of  roads  in  the  neighborhood 
of  the  contracter's  plant.  The  GGSS  van  will  be  driven  in  a  test  pattern  similar 
to  those  use  to  test  an  inertial  postltloning  system.  After  the  vehicle  has 
completed  a  number  of  repeat  tracks,  a  study  will  be  made  to  determine  instrument 
parameters  and  sensitivities.  Noise  sources  are  then  Isolated  and  corrected  or 
compensated  for  In  the  software. 


Similar  repeat  tracks  are  then  flown  with  the  van  in  the  aircraft.  In  each 
case,  parameters  such  as  vehicle  velocity,  track  control  mechanisms  (autopilot, 
cruise  control,  and  manual)  and  road  surface  or  atmospheric  turbulence  will  be 
varied  to  test  for  Instrument  sensitivity  to  a  particular  parameter. 

While  the  statistical  characteristics  of  the  gradient  field  will  show  up  in 
this  test,  it  is  not  critical  that  they  be  understood  exactly  at  this  point  in 
testing.  The  test  method  here  is  repeatability.  Therefore,  it  is  only  important 
that  the  gradient  field  have  no  temporal  variations,  so  that  the  difference 
between  each  traverse  can  be  judged  as  an  error. 


Phase  II  Testing 


During  Phase  II  testing,  each  vehicle  will  "survey"  a  predetermined  set  of 
tracks  and  be  required  to  produce  estimates  of  gravity  at  points  along  the 
track  in  the  land  case,  and  at  selected  points  beneath  the  flight  grid  for  the 
airborne  test.  During  this  phase  of  testing,  the  system's  performance  will  be 
evaluated  under  near  "optimum"  conditions,  subject  only  to  changes  In  the 
strength  of  the  gravity  field.  This  will  allow  an  understanding  of  the 
gradiometer 's  ability  to  recover  the  gravity  f^cJd  under  optimal  "operational" 
conditions.  Both  a  relatively  smooth  and  a  relatively  rough  gravimetic  field 
are  found  in  the  Phase  II  test  area  (Figure  9). 
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Phase  II  Land  testing  will  verify  the  GGSS's  ability  to  "transfer"  gravity 
and  deflection  of  the  vertical  between  a  point  with  a  known  gravity  vector  and  a 
point  with  unknown  values.  A  series  of  increasingly  longer  loops  will  be 
traversed  starting  at  the  Cllnton-Sherman  Airfield  in  Burns  Flat,  Oklahoma 
(Figure  10).  On  each  of  these  loops  a  number  of  "truth  data"  points  will 
be  presurveyed.  Comparisons  will  be  made  at  a  selected  number  of  these  points. 

The  mismatch  between  "truth"  and  the  measured  value  will  be  a  direct  indication 
of  the  error  in  the  GGSS's  recovery  of  the  gravity  vector. 

Phase  II  Airborne  testing  will  verify  the  capability  of  the  GGSS  to  provide 
survey  quality  gravity  data  from  an  airborne  platform.  The  GGSS  van  will  be 
loaded  onboard  the  C-130  aircraft  and  then  fly  a  grid  pattern  like  that  of 
figure  11.  An  area  of  approximately  three  hundred  kilometers  square  will  be 
covered  with  orthogonal  tracks  at  five  kilometer  track  intervals  by  the  end  of 
the  test  program. 

A  day’s  survey  will  be  controlled  by  the  operator  through  the  airborne 
computer.  Flight  planning  information  will  be  stored  by  the  operator  prior  to 
each  day's  mission  in  the  form  of  "waypoints".  Once  the  aircraft  is  in  the 
survey  area,  the  GGSS  computer  will  fly  the  aircraft  through  an  autopilot 
interface.  A  set  of  cockpit  Informational  displays  and  Instruments,  and  control 
swltchs  will  allow  the  pilot  to  engage,  disengage  and  monitor  the  progress  of 
the  survey. 

A  typical  flight  will  last  six  to  eight  hours  at  a  constant  altitude  near 
six  hundred  meters  above  ground  level  and  at  a  constant  ground  speed  of  about 
four  hundred  fifty  kilometers  per  hour.  Operation  of  the  GGSS  requires  that  at 
least  three  GPS  satellites  are  visible  while  the  survey  is  in  progress.  Because 
GPS  will  only  be  partially  operational  during  this  test  program,  GGSS  testing 
can  only  be  conducted  only  in  a  satellite  availability  "window".  With  the  present 
constellation  of  satellites,  approximately  six  hours  of  coverage  will  be  available 
each  day.  The  number  of  satellites  and  the  times  of  day  they  will  be  visible 
during  October  of  86  is  shown  in  figure  12.  The  general  shape  of  this  bar 
chart  does  not  change  from  day  to  day.  £ach  month,  the  "window"  moves  about  two 
hours  earlier  in  time. 

It  is  the  philosophy  of  this  test  program  to  attempt  to  measure  the  Instruments 
potential  for  providing  needed  gravity  survey  capability.  Phase  11  testing  is 
specifically  designed  with  as  few  obscuring  variables  as  is  possible.  Many 
factors  can  and  will  degrade  system  performance.  The  "noisy"  aircraft  vibration 
and  turbulence  environment,  poor  aircraft  control  parameters:  altitude  hold 
times,  poor  tracking  or  airspeed  control  or  uncoordinated  flight,  and  suboptlmal 
data  processing  are  all  examples  of  possible  degrading  factors.  Others  will  be 
found  in  the  course  of  testing.  Like  other  instrument  tests,  one  tries  to 
control  as  many  of  the  known  parameters  as  possible  until  they  are  characterized 
or  their  influence  mitigated.  Compensation  for  those  that  cannot  be  controlled 
will  be  accomplished  post  mission  using  signal  processing. 

For  the  land  vehicle,  an  attempt  will  be  made  to  maintain  an  optimum  constant 
speed  on  relatively  "good"  roads,  free  of  interfering  traffic.  Events  such  as 
the  passing  of  a  large  truck  will  be  logged  for  analysis  and,  if  needed,  later 
compensation. 


In  the  airborne  environment,  every  attempt  will  be  made  to  minimize 
Interference  from  other  air  traffic.  The  test  will  be  flovm  under  near  perfect 
aviation  conditions.  Altitude  and  airspeed  will  be  maintained  near  constant, 
and  heading  will  remain  constant  for  long  periods  of  time.  The  aircraft's 
airspeed  will  be  maintained  so  as  to  minimize  vibration  and  turbulence  induced 
acceleration  errors  found  in  Phase  I  testing  .  Flights  will  probably  be  flown 
during  or  at  least  with  prior  knowledge  of  optimal  GPS  satellite  geometry.  A 
grid  of  evenly  spaced  tracks  will  be  flown  with  every  attempt  to  maintain  the 
uniformity  of  the  grid.  Attempts  have  been  made  to  provide  a  flight  environment 
which  minimizes  weather  induced  errors  and  delays,  air  traffic  conflicts  and 
deviation  from  the  constant  set  of  test  conditions.  Therefore,  the  test 
should  conclude  with  an  accurate  forecast  of  the  potential  value  of  the  GGSS  to 
a  survey  organization,  like  DMA. 

Laboratory  testing  will  begin  in  January  of  1986  with  Phase  I  and  II  to 
follow.  Initial  testing  is  expected  to  be  completed  in  early  1987.  A  summary 
of  the  test  program  is  provided  in  Table  2.  Follow-on  tests  are  anticipated 
and  will  evaluate  the  best  way  to  deal  with  operational  problems  such  as  variable 
altitudes,  airspeeds,  topographic  relief  patterns,  weather,  and  aircraft 
capabilities . 

The  development  of  the  GGSS  is  being  managed  by  the  Earth  Sciences  Division 
of  the  Air  Force  Geophysics  Laboratory  (AFGL),  Hanscom  AFB,  Massachusetts  under 
funding  from  the  Defense  Mapping  Agency  (DMA),  Washington  D.C..  The  GGSS  is  a 
product  of  Bell  Aerospace  Textron,  Niagara  Falls,  New  York. 
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Figure  1  -  This  figure  shows  two  accelerometers  located 
side  by  side  and  with  sensitive  axes  aligned  in  the  s 
direction.  Each  accelerometer  Is  measuring  the  total 
force  applied  to  it.  Both  accelerometers  are  subject  to 
an  inertial  acceleration  a  and  gravitational  force  which 
differs  byiig  over  Che  distance  seperaCing  them^x.  The 
difference  in  applied  force  divide  by  the  distance 
separating  the  accelerometers  is  Che  measurement  of  the  t 
component  of  Che  gravity  gradient  in  the  x  direction.  In 
ocher  words,  it  is  the  difference  in  the  z  component  of 
gravity  over  Che  distance  Ax. 
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Figure  2  -  The  gravity  gradient  produced  by  a  mass 
located  at  the  given  distance  from  a  gradioineter 
is  plotted  here.  Examples  of  the  gradient  produced 
by  masses  of  Interest  in  the  testing  of  a  gravity 
grad tome ter  are  shown. 
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Figure  3  -  Bell  Model  VIIB  Accelerometer 


Figure  5-The  GGI  is  the  heart  of  the  GGSS  system. 

It  contains  four  of  the  BELL  Model  VII  accelerometers, 
mounted  tangentially  to  the  edge  of  a  rotating  disk, 
at  90  degree  intervals  .  The  sum  of  the  outputs  of  a 
set  of  accelerometers  gives  the  average  gradient  across 
the  distance  of  separation.  A  set  of  three  GGIs  are 
mounted  on  the  gryo  stabilized  platform  shown. 


FIGURE  6  - A  SIMPLIFIED  FUMCTIOIIAL  BLOCK  AND  SIGNAL  FLOW  DIAGRAM  FOR  THE  GGSS 
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STAGE  2  AIRBORNE  PROCESSING 
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Figure  7-Block  Diagram  of  Che  Stages  of  Data  Processing 
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Figure  9  -  This  aep  shows  the  Fhsse  II  GG88  Test  Ares. 
The  outline  shows  the  300  by  300  kllosKter  sres  lAIbh 
will  be  flown  by  the  GG8S  ss  a  test  of  surwey  capability 


AIRBORNE  GRAVITY  GRADIOMETRY 


Figure  11  -  Schematic  of  Phase  II  Aircraft  Testing. 
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Figure  12  -  The  Global  Posicioning  Syacem  (GPS)  will  be  the  prinary 
navigation  tool  used  during  this  test  prograa.  Four  visible  satellites 
are  required  to  obtain  a  three  diisensional  position  without  other  aids 
such  as  an  atonic  clock  or  barometric  altineter.  GPS  will  provide 
worldwide  navigation  capability  in  the  near  future.  The  present 
constellation  of  satellites  provides  four  satellites  visibility  only 
a  United  nunber  of  hours  per  day.  The  projection  for  an  average  day 
in  October  1986  is  given  above.  Tlie  height  of  the  bar  indicates  the 
number  of  satellites  thst  will  be  visible  at  that  hour  of  the  day  in 
the  Texas /Oklahoma  Area. 
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ANTICIPATED 


GRAVITY  DISTURBANCE  RECOVERY 
(WAVELENGTHS  <  500  KM) 

DEFLECTION  OF  THE  VERTICAL 

POSITION  ERROR 

PLATFORM  TILT 

AZIMUTH 

SYSTEM  NOISE  (maximum  allowable) 


.9  fflgal 

.83  mgal 

.18  arcsec 

.16  arcsec 

100  meter 

15  meters 

26  arcsec 

15  arcsec 

66  arcmln 

34  arcmln 

10“3  +  190  (E/R/S) 


TABLE  1  -  AIRBORNE  GGSS  ACCURACY  REQUIREMENTS 


Test 


Method 


Purpose 


LABORATORY 

PHASE  I 

LAND  & 
AIRBORNE 

PHASE  II 

LAND 

AIRBORNE 


LABORATORY  METHODS 

Study  lastrument  outputs 
produced  by  various  Inputs 
under  static  and  dynamic  but 
controlled  conditions. 

REPEATABILITY  TEST 

Comparison  of  the  ability 
of  the  Instrument  to  give 
consistent  results  when  the 
Instrument  traverses  the 
same  test  track  many  times. 

COMPARISON  TO  "TRUTH"  DATA 


Comparison  to  "truth"  data 
by  along  track  Integration 
between  a  "truth"  data  point 
and  a  test  point. 

Comparison  to  "truth”  data 
by  gathering  data  In  a  large 
area  and  reducing  this  data 
to  gravity  values  at  selected 
point  where  "truth"  data  Is 
aval lb le . 


Calibration  and  compensation 

Understand  Instrument  outputs 
and  their  sensitivities 


Verify  Instrument  noise 
and  sensitivity 

Characterize  the  statistical 
nature  of  the  gradient  field 


Verify  the  system  Is  capable 
of  providing  survey 
quality  data 


Table  2  -  GGSS  Testing  in  Summary 


3 


AIRBORNE  GRAVITY  GRADIENT  SURVEY  DATA  REDUCTION 

W.  J.  Hutcheson 
Bell  Aerospace  Textron 
P.  0.  Box  One 
Buffalo,  NY  142no 


ABSTRACT 

Algorithms  which  have  been  developed  at  Bell  Aerospace  Textron  for  the 
GGSS  Phase  II  test  plan  data  reduction  requirements,  will  be  described  in 
detail.  The  data  reduction  stages  include  track  crossing  adjustment,  optimal 
spatial  integration,  gridding.  Interpolation  and  smoothing,  low  frequency 
adjustment  using  astrc^eodetic  data  and  downward  continuation.  Simulation 
results  obtained  using  a  synthetic  field  will  be  {presented  to  demonstrate 
the  algorithm  performance. 
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AVERAGE  POSITIVE  ERROR  OF  CENTRAL  STRIP  PASSING  THROUGH  FIELD 


^EOyENCY  ERROR  ESTIMATED  USING  TRUTH  DATA  FIELD  LOW  FREQUENCY  DATA  ESTIMATED  USING  9  ASTROS 


HORIZONTAL  DISTURBANCE  AVERAGE  ERROR  OVER  GRID  SQUARE 


Bell  Aerospace 


PAPER  TITLE:  AIRBORNE  GRAVITY  GRADIENT  SURVEY  DATA  REDUCTION 

SPEAKERS  NAME:  John  Hutcheson 


Questions  and  Comments: 

•  Julian  Center:  What  statistical  model  do  you  use  for  least  squares  collocation 
to  incorporate  tie  points? 


SPEAKERS  RESPONSE:  The  STAG  (or  reciprocal  distance)  7-tenn  gravity  model 
is  used  in  the  LSC  algorithm  to  upward  continue  the  astro  tie  point  to  a 
node  on  the  derived  map. 


ISOMORPHIC  GEODETIC  AND  ELECTRICAL  NETWORKS:  AN  APPLICATION  TO 


THE  ANALYSIS  OF  AIRBORNE  GRAVITY  GRADIOMETER  SURVEY  DATA 


RRH:  GRAVITY  GRADIOMETER  SURVEY  ANALYSIS 


D.  H.  Eckhardt 

Air  Force  Geophysics  Laboratory 
Earth  Sciences  Division 
Hanscom  AFB,  MA  01731-5000 

ABSTRACT 

Late  In  1986,  the  Bell  Aerospace /Textron  Gravity  Gradlometer  Survey  System 
(GGSS)  will  be  tested  by  the  Air  Force  Geophysics  Laboratory  In  an  airborne  survey 
of  a  300  km  x  300  km  region  of  Oklahoma  and  Texas*  The  survey  pattern  will  be  a 
grid  with  a  5  km  separation  between  adjacent  tracks,  north-south  and  east-west. 

One  way  to  process  the  GGSS  survey  data  Is  to  analyze  an  electrical  network  that  is 
Isomorphic  to  the  survey  network.  The  Integrated  gradients  between  the  survey 
crossing  nodes  correspond  to  the  applied  voltages  between  the  nodes  of  the  electri¬ 
cal  network;  the  gradient  variances  correspond  to  the  inter-nodal  resistances;  the 
elements  of  the  adjusted  gravity  vector  correspond  to  the  nodal  voltages;  and  the 
solution  variances  correspond  to  the  resistances  to  ground.  The  survey  error 
analysis  Is  performed  by  calculating  the  resistance  to  ground  of  the  electrical 
network;  a  technique  for  making  the  calculations  in  large  networks  is  explored  in 
detail.  For  a  sample  survey  scenarios  with  one  ground  truth  control  point  near  each 
corner  of  the  survey  square  and  with  realistic  values  for  the  survey  parameters, 
the  gravity  disturbance  standard  deviation  Is  less  than  0'.'25  at  all  nodes.  With  no 
'ground  truth,  but  with  a  gravimeter  on  the  aircraft  that  can  independently  determine 
gravity  to  10  mGal  at  all  nodes,  the  adjusted  gravity  disturbance  standard  devia¬ 
tion  Is  less  than  1  mGal. 


FI(J.  6.  This  Infinite  ladder  of  resistors  Is  a  pantopod  network 


ladder  of  resistors 


FIG.  10.  TMs  Infinite  lattice  of  resistors  Is  a  nonopode  network 


I!,  'tip  (0.10  mf.al  contour  Int^rvh.!)  of  gravity  disturbance  standarl 
l»'vIitlon  In  test  survey  area;  at  the  four  tie  points  near  the 
'•arners  the  standard  deviations  are  0.46  aCal  a  priori,  and 
0.'«2  mOal  a  posteriori. 
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STAGE  II  PROCESSING  OF  AIRBORNE  GRAVITY  GRADIOMETER 
DATA  USING  FREQUENCY  DOMAIN  TECHNIQUES 
A.  A.  Vasfiiliou 
The  University  of  Calgary 
Division  of  Surveying  Engineering 
2500  University  Drive  N.W. 

Calgary,  Alberta,  T2N  1N4 
Canada 


ABSTRACT 

A  set  of  frequency  response  functions  between  the  first-order  and  the 
second-order  gravity  gradients  Is  developed  In  this  paper  using  flat-earth 
approximation.  In  this  way  the  spectrum  of  each  first-order  gradient  is 
^related  to  the  spectra  of  one  or  more  second-order  gravity  gradients.  Assuming 
noise-free  second-order  gradient  measurements,  these  equations  are  transformed 
back  Into  the  space  domain  as  Integral  equations  relating  each  first-order 
gradient  to  all  of  its  second-order  gradients  (e.g.  Tr  to  Txz,  Tyz,  Tzz). 

The  frequency  domain  relations  are  used  to  estimate  the  first-order  gradients 
employing  FFT .  A  set  of  simulated  data  Is  used  to  test  these  relationships. 
Results  from  these  tests  and  computer  time  requirements  are  finally  discussed 


In  this  paper. 


STAGE  II  PROCESSING  OF  AIRBORNE  GRAVITY 


GRADIOMETER  DATA  USING  FREQUENCY  DOMAIN  TECHNIQUES 


Anthony  A.  Vassillou 


The  University  of  Calgary 
Division  of  Surveying  Engineering 
2500  University  Drive  N.W, 
Calgary,  Alberta  T2N  1N4 
CANADA 


the  14th  Annual  Gravity  Gradlometer  Conference 


Presented  at 


Colorado  Springs,  Colorado 
February  11-12,  1986 


OBJECTIVES : 


1.  Develop  a  fast  computational  method  in  the  frequency  domain  for  the 
estimation  of  the  first-order  anomalous  gravity  gradients  from 
gravity  gradlometer  data.  Take  advantage  of  the  efficiency  of  FFT. 


2.  Use  as  many  second-order  gravity  gradients  as  possible  for  the 
estimation  of  any  single  first-order  gradient.  Take  the 
gradlometer  self-noise  explicitly  Into  account. 


DATA  USED: 


A  set  of  about  250  000  simulated  airborne  gradient  measurements  on 
a  2.31  X  2.44  km  grid,  over  a  472  (north)  x  496  (east)  km  area.  Is 


used . 
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2.  INPUT  -  OUTPUT  FILTERING  EQUATIONS 


2.1  Single  Input  -  Single  Output  Equations 


f (x,y) 


h(x,y) 


g(x>7) 


G(u,v)  »  H(u,v)  F(u,v) 


S  (u,v)  *  1h(u,v)|  S,,(u,v) 

gg  If 


Sj^(u,v)  *  H(u,v)  Sj^(u,v) 


For  extraneous  noise  n(x,y)  present 


S  (u,v)  *  1h(u,v)|  {S,,(u,v)  +  S  (u,v)} 

gg  If 


*  H(u,v)  {S,,(u,v)  +  S  (u,v)} 
1 1  nn 


S,  (u,v) 
fg 


(5) 
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2.2  Multiple  Input  -  Single  Output  Filtering  Equations 

For  simplicity  reasons •  discuss  only  two  input  -  single  output 
systems 


f ^(x,y) 


h2(x,y) 


E 


>  g(x,y) 


•  The  two  Inputs  and  output  related  in  the  frequency  domain  by 

G(u,v)  »  Hj(u,v)  F^(u,v)  ♦  H^Cu.v)  F^fu.v)  (6) 

•  For  partially  correlated  Inputs  fj(x,y),  f2(x.y),  the  transfer 
functions  are  given  by 


S,  ,  (u,v)  S,  „(u,v)  -  S,  ,  (u,v)  S,  „(u,v) 

.  ^2’^2  ^'8  ^1*^2  ^2*8 

H.(u,v)  =  - 1 - s -  (7) 

S,  ,  (u,v)  S,  ,  (u,v)  -  |S^  ,  (u,v)! 


'fj.fl 


^2’^2 


f,.f2 


S  (u,v)  S-  „(u,v)  -  S,  ,  (u,v)  (u,v) 

*2*8  *2’*!  *r8 

H2(u.v) - ^ ^ -  (8) 


fj.fj 


(u,v)  S,  ,  (u,v)  -  Is,  ,  (u,v)!^ 


'f2.f2 


fj,f2 


For  fully  correlated  Inputs  ^^*^2’  transfer  functions  ^^*^^2 
the  form 


H^(u,v) 


S  _^u,v) 

^'8 


S-  ,  (u,v)  ♦  S,  (u,v) 

tl.t,  t2'^2 


(9) 


H^Cu.v) 


S  ,  (u,v)  +5,  f  (u,v) 
1’  1  2*  2 


(10) 
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For  input  noise  uncorrelated  to  both  fully  correlated  input  signals, 
the  transfer  functions  are  modified 


Hj(u,v) 


H^Cu.v) 


S,  (u.v) 

fi.g 


S  ,  (u.v)  +  S  (u.v)  +  S  „(u,v) 

fj.fj  f2»t2 


S-  _(u.v) 

^2 


S  ,  (u.v)  +  S  -  (u.v)  +  S  „(u.v) 
fj^.fj  ^2*^2  n.n 


(11) 


(12) 


Similar  expressions  to  (9),  (10),  (11),  (12)  can  be  derived  for  a 
case  of  more  than  two  fully  correlated  input  signals. 


Multi-input  -  Single  output  filtering  equations  very  useful  for 

airborne  gravity  gradlometer  where  all  input  signals 

T  ,  T  ,  T  are  fully  correlated, 
yy  zz  yz 
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3.  APPLICTION  OF  MULTIPLE  INPUT  -  SINGLE  OUTPUT  SYSTEMS  TO  THE 

ESTIMATION  OF  FIRST-ORDER  GRADIENTS  FROM  SECOND-ORDER  GRADIENT  DATA 

Statement  of  the  problem:  Estimate  T  from  T  ,  T  ,  T  ,  T  ,  T 

- -  z  xz  yz  zz  XX  yy 

Flat  earth  approximation  Is  assumed. 


3a.  Estimation  of  T  from  T  gradient  measurements 

Z  Z  Z 


•  To  estimate  T  (x,y)  from  measurements  of  T  (x,y,  z»h) ,  is 
z  zz 

necessary  to  derive  the  transfer  function  between  them 


H(u,v) 


T  ,T 


(u,v)+S  (u,v) 

n  an 


where 


(u' 


Sx^n^e  _(u,v) 


T,T' 


Ibx^q^e  _(u,v)+S  (u,v) 

T  ,T  n ,n 


(13) 

(14) 


For  simplicity  reasons  assume  noise-free  measurements 
2xhq 

H(u,v)  =  -  - -  (15) 

q 

T^(x,y)  cannot  be  computed  at  the  (0,0)  point.  The  mean  of  T^  has  to 
be  evaluated  from  other  sources. 


•  Equations  (13),  (14)  incorporate  plane  integration  and  downward 
continuation. 

Pure  plane  integration  corresponds  to  Stokes  formula 
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3b.  Estimation  of  T  from  T  -  T 

z  xz  yz 


Necessary  to  determine  the  transfer  functions  between  T  -  T 

’  xz  yz 

and  T  (noise-free  measurements  assumed) 
z 


Hj(u,v) 


H2(u,v) 


P{T^(x,y)}  -  H^(u,v)  P{T^^(x,y  z-h) }  +  H^Cu.v)  F{T^|x,y,  z=h) } 


where  F  denotes  the  Fourier  transform. 


*  Pure  plane  integration  transformed  in  the  space  domain  yields 


-  -  H  - - ^2- - '’’'l'*’'!-  2V  ,  ^  ,2^, - ^ 

[(x-Xj)  ♦(y-y^)  y  [x-x^)  ♦(y-y^)  ] 


where 


(x-x^ 


2  2  is  ’  COSO  -  „  ^  . 

((x-Xj)^+(y-yj)^]’*  [{x-Xj)S(y-y^)n^ 


(y-Fi) 


Equation  (20)  Is  Important,  providing  directly  in  terms  of 

A  similar  equation  can  be  derived  for  the  anomalous  potential 


T  (x,,y,)  slna 

[(x-Xj)^(y-y.)^P  ^  J  J  [(x-x^)^(y-y^)-]^ 


According  to  (22)  ,  geoidal  undulations  can  be  computed  in  a  local  area 


from  gridded  2-D  deflections  of  the  vertical  (ideal  case). 
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3c.  Esclmation  of  T  from  T  -  T  -  T 
_ z _ xz  yz  zz 


I 


•  Determine  transfer 
free  measurements 


functions  between  T  -T  -T  and  T  for  i.olse 

xz  yz  zz  z 


H^(u.v) 

H2U,v) 

H3(u,v) 


2xhq 


1 


2xhq 
u  ♦  e  ^ 

2Trq^ 

^2xhq 

4xq 


(23) 

(24) 

(25) 


F{T  (x,y)}  *  H  (u,v)P{T  (x  • .z-h)>+H*(u,v)F{T  (x,y,z»h)}+H  (u,v) 

z  X  xz  ^ 

P(T^^(x,y,z-h)}  (26) 

•  Pure  plane  Integration  transformed  In  the  space  domain  results  In 


T^(x.y)  - 


\_ 

4ti 


T^^(x,,y,)  slno 
xz  1  1 


'‘’l'*’'!'  47 


r 


T  (x,,y,)  COSO 
yz  1  -^1 


)  } 


I(x-Xj)‘4(y.yj)‘) 


1_ 


.OD 


I(x-Xj)%(y-yj)‘l 


IS 


T  related  directly  to  T  ,  T  ,  T  by  an  integral  formula, 
z  ^  xz  yz  zz  * 


•  A  similar  formula  can  be  developed  for  the  anomalous  potential  T 


(27) 


T  (x.y) 


00  cc 


l(x-xj)^*(y-yj)^]** 


dxjdyj- 


4r, 


|(x-Xj)^.(y-y|)^l** 


dxjdy, 


411 


(x,.y,) 


!(x-x,)‘.(y-yj)‘| 


IS  "’‘l'*''! 


•  Formula  (28)  provides  the  anomalous  potential  (g^oldal  undulations)  In 
terms  of  all  of  its  first-order  gradients.  For  local  gravity  approxi¬ 
mation  though,  this  formula  is  not  applicable  due  to  the  unavaila¬ 
bility  of  T  ,T  (deflection  data)  in  2-D  grids.  Nevertheless  equation 
(28)  is  very  useful  because  all  second-order  gradients  T  ,  T  ,  T 
are  available  from  a  gravity  gradiometer  system.  ^ 


(28) 
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Notes: 


1. 


In  sections  3a.  3b.  3c. 
substituted  by  the  sum  - 
form 


the  second-order  gradient  T  can  be 

(T  +1  ).  So  equation  (f5)  takes  the 

XX  yy  ^ 


1  2xhq 

F{T^(x.y)}  -  — P{T^^(x,y,z-h)  +  (x.y ,z-h) } 


(29) 


2.  All  the  multiple  Input  -  single  output  equations  are  equivalent  to 
the  ones  derived  from  a  multidimensional  Wiener  filtering  with 
linearly  correlated  Inputs.  Wiener  filtering  (in  a  planar 
approximation)  can  be  formulated  as  the  2-D  Fourier  transform  of  a 
least-squares  collocation  estimate.  However,  Wiener-filtering  is 
much  faster  computationally  to  least-squares  collocation.  For 
Instance  assuming  a  2-D  grid  with  ^  points  In  each  direction  and 
three  gradient  measurements  per  point,  least-squares  collocation 
has  to  Invert  a  22500  x  22500  matrix.  This  takes  more  than  15  CPU 
hours  on  the  Supercomputer  Cyber  205.  On  the  other  hand 
application  of  multi-dimensional  Wiener  filtering  or  multiple 
input-output  filtering  equations,  leads  to  an  inversion  of  a  3  *  3 
matrix  for  each  pair  of  frequencies  (u.v).  For  linearly  related 
Inputs  this  inversion  reduces  dovm  to  an  inversion  of  a  1  «  1 
matrix. 


3.  So  far  the  effect  of  topography  on  the  gradients  has  been 

neglected.  However  it  can  be  explicitly  computed  from  the  spectra 


XX 


yy 


zz 


xy 


xz 


yz 


S.^nf-1  h"(x.y)} 

(30) 

n*l 

n!  q 

2xGe"^’^‘’*o  1 

—  2irjr_  h"(x,y)} 

(31) 

n*l 

n!  q 

2:TGe"^’^‘*^o  1 

^213)1  F{p(x.y)  h"(x.y)} 

(32) 

n«=l 

n! 

2nGe"^’'‘*^o  Z 

“’"i’"-.  2*“''  F{(,(x.y)  h"(,.v)) 

(33) 

n“l 

n!  q 

2irGe"^^‘**o  E 

- j  2iru  F{p(x.y)  h”(x,y)} 

(34) 

n*l 

n! 

2irGe'^’'‘’*o  E 

- 3 -  j  2irv  F{p(x,y)  h  (x,y)} 

(35) 

n«l 

n? 
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Note  3.  (continued) 


where  in  equations  (30)  to  (35)  F{T^j }  represent  the  spectrum  of  the 


effect  of  topography  on  the  second-order  gradient  T 
z  is  the  flying  altitude 


IJ’ 


h(x,y)  is  the  grldded  height 

p(x,y)  is  the  gridded  density  (it  can  be  taken  as  constant 

2.67  gr/mh 
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4.  DATA,  TESTS  AND  RESULTS 
Aa  Simulated  data 

•  Airborne  and  earth's  surface  first-order  and  second-order  gravity 
gradients  are  simulated  on  a  2.31  km  (north)  *  2.44  km  (east) 
grid.  For  the  simulation  9.26  km  (north)  *  9.74  km  (east) 

(i.e.  5'  *  10’)  gravity  anomaly  data  are  used.  The  gravity 
anomaly  data  are  given  on  a  grid  In  Northern  Saskatchewan,  Canada. 
Extent  of  the  grid  Is  472  *  496  km. 


•  Flat-earth  approximation  Is  used  under  the  transformation 
dx  *  R  cos^  dX 


R  -  6371  km 


dy  *  Rd^ 


•  Long-wavelengths  are  eliminated  by  subtracting  the  spherical 
harmonic  expansion  up  to  degree  and  order  36  of  Rapp  1978  model, 
from  the  original  free-air  gravity  anomalies. 


•  The  modelling  of  the  gravity  data  Is  made  by  a  two-layer  point 
mass  model  (depths  at  15.5  km  and  5  km  respectively).  For  high 
frequency  information  a  third  layer  Is  added  burled  at  depth  of 
1  km,  with  whlte-nolse  point  mass  distribution  on  it. 


The  second-order  gradients  are  simulated  at  a  flying  altitude  of 
600  m.  They  were  corrupted  by  gradlometer  noise  (following  the 
model  published  by  White  (1980)). 
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4b.  Results  from  the  estimation  of  T  from  T  ,  T  ,  T 

Z  X2  yz  22 


The  mean  of  T  over  the  whole  area  is  assumed  known  from  other 
z 

sources.  Results  are  given  on  the  Inner  area  bounded  45  km  from 
the  north  and  the  south  borders  and  48  km  Inside  the  east  and  west 
borders. 


472  km 


1st  Case 


48  k 


m 


Is  estimated  directly  below  measurement  grid  points 


Measurements 


T 

zz 

T  -  T 
xz  yz 

T  -  T  -  T 
xz  yz  zz 

RMS  error  of  T 

z 

0. 70  mgals 

0.46  mgals 

0.41  mgals 

CPU  time 

(on  a  Honeywell 
Multlc  DPS  68 
computer) 

545  CPU 
seconds 

560  CPU  seconds 

578  CPU  seconds 
(124848  measurements) 
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2nd  case  :  The  gradlometer  self-noise  is  taken  into  account. 

The  PSD  of  the  anomalous  potential  T  is  modelled  by 


where  A, 
Saskatcn 


has  been  determined  from  data  of  all  the 
evan  province. 


Measurements 


T 

zz 

T  -  T 
xz  yz 

T  -T  -T 
xz  yz  zz 

RMS 

error 

of  T 

2 

0.63  mgals 

0.44  mgals 

0.37  mgals 

CPU 

time 

622  CPU  seconds 

640  CPU  seconds 

660  CPU  seconds 

1st 

case 

:  The  gradlometer  self-noise 

isn’t  taken  into 

account . 

T  is 
z 

estimated  by  FFT  interpolation  at  all 

the  points  on 

1.155 

X  1.22  km  (north-east)  grid. 

Measurements 

T 

zz 

T  -T 
xz  yz 

T  -T  -T 
xz  yz  zz 

RMS 

error 

of  T 

z 

0.50  mgals 

0.38  mgals 

0.32  mgals 

(estimation  at 
499  392  points) 
1  680  CPU  seconds 


CPU  time 


lit 


BLOCK  DXAORAH  OF  T2  QBAOICNTS 


PLOT  NO.  1 
A2IM  •  -85. 0 


DATE  18<<19/85 
ELEV  •  8S.0 


time  19. 
DZ8T  > 


11  .5 
10000 


BLOCK  DIAGRAM  OF  T2  ER.  FROM  TZ2 
PLOT  NO.  1  date  la/eS^BS 

AZIM  «  -8S.0  ELEV  «  8S,0 


time  19.09.S 
DIST  •  10000 


TYZ 

time:  19.16.3 
OIST  »  10000 
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CONCLOSIONS 


A  fast  computational  method,  employing  FFT,  Is  developed  for  the 
estimation  of  first-order  gradients  from  second-order  gradients. 
The  method  Is  based  on  the  application  of  the  multiple  Input  - 
single  output  filtering  equations. 


The  method  handles  all  possible  combinations  of  gradlometer  data. 

For  the  estimation  of  T  uses  T  ,  T  -T  ,  T  -T  -T 

z  zz  xz  yz  xz  yz  zz 

Furthermore  takes  the  noise  explicitly  Into  account. 


Results  from  test  runs  with  simulated  data  show  that  T  can  be 

z 

estimated,  directly  below  the  gradlometer  points  with  an  RMS  error 
of  less  than  0.70  mgals.  The  points  at  the  middle  of  the  gradlo¬ 
meter  points,  can  be  Interpolated  with  an  RMS  error  of  less  than 
0.5  mgals. 


il 


KARHUNEN-LOEVE  GRAVITY  GRADIOMETER  DATA  PROCESSING 

S.  Bose 

Applied  Science  Analytics,  Inc. 

7041  Owensfflouth  Ave. 

Suite  206 

Canoga  Park,  CA  91303 


ABSTRACT 

The  Karhunen-Loeve  data  processing  technique  for  gravity  gradlometry  Is 
reviewed.  The  method  Is  based  on  the  White  Noise  Layer  (WNL)  ^jodel  which  replaces 
the  unknown  mass  distribution  below  the  survey  region  as  multiple  two-dimensional 
white  noise  layers  representing  the  vertical  derivative  of  the  disturbance  potential 
to  any  order.  Such  a  gravity  signal  model  is  derived  from  the  physical  theory  of 
geodesy  and  Is  particularly  suited  for  modeling  high  frequency  phenomena.  Of 
particular  Interest  Is  that  such  a  model  results  in  a  nonstationary  nonlsotroplc 
representation  of  the  disburbance  potential.  A  two-dimensional  signal  processing 
algorithm  to  process  all  the  gradlometer  data  simultaneously  is  presented.  The 
estimation  algorthm  can  handle  multiple  layers  of  different  inline  or  crossline 
data  given  In  two  dimensional  grids  at  the  same  or  different  altitudes  on  or  above 
the  surface  of  the  earth.  Different  grid  patterns  of  the  same  Inline  or  crossline 
measurements  at  the  same  or  different  altitudes  are  also  accommodated  by  the 
algorithm  without  resorting  to  any  pre-processor  averaging  techniques.  The 
method  Is  such  that  at  any  given  spatial  point  the  gravity  field's  correlation  in 


any  direction  Is  not  ignored  and  the  estimation  algorithm  does  not  enforce  any 
unnecessary  limitation  of  casuality  on  the  data  Inasmuch  as  no  one-dimensional 
scanning  is  performed.  The  problem  of  simultaneous  smoothing  of  all  the  gradio- 
meter  measurements  from  all  survey  traverses  in  the  region  is  solved  by  represent¬ 
ing  the  disturbance  potential  solution  as  a  Karhunen-Loeve  (KL)  expansion. 
Estimating  the  gravity  field  or  any  of  its  derivatives  simply  reduces  to  estimat¬ 
ing  the  KL  coefficients  and  performing  the  appropriate  transformations  involving 
sine  and  cosine  functions.  Neither  the  estimation  of  the  KL  coefflcents  nor  the 
implementation  of  these  transformations  requires  any  matrix  inversions.  Dotmward 
continuation  as  well  as  Interpolation  of  estimates  between  grid  measurements  are 
performed  automatically.  Thus  large  amounts  of  moving  base  gravity  gradiometry 
data  can  be  handled  by  this  technique  in  a  computationally  efficient  manner. 
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APPLIED  SCIENCE  ANALYTICS,  INC. 


DEFINITION  OF  TRANSFORMATION  MATRICES 


APPLIED  SCIENCE  ANALYTICS,  INC 


UTILIZATION  OF  SPECIAL  PROPERTIES  OF  A  TOEPLITZ  CIRCULANT  MATRIX 
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ESTIMATION  ALGORITHM  REQUIRES  NO  MATRIX  INVERSIONS 
MEASUREMENT  DATA  MUST  BE  IN  PLANAR  GRIDDED  FORM 


PAPER  TITLE: 


KARHUNEN-LOEVE  GRAVITY  GRADIOMETER  DATA  PROCESSING 


SPEAKERS  NAHE:  Sam  Bose 
Questions  and  Comments: 

John  V.  Breakwell:  Aren't  you,  in  effect,  assuming  that  your  finite  gravitational 
ground  pattern  Is  repeated  (to  Infinity  In  all  4  directions)? 

SPEAKERS  RESPONSE :  No 

Julian  Center:  1.  Why  are  there  only  sine  terms  In  the  K-H  expansion  for 
the  potential? 

2.  How  does  your  metnod  differ  from  the  Pour ler-Trans form/  Wiener  Smoother 
Method? 

SPEAKERS  RESPONSE:  1.  Only  sine  terms  in  representation  of  T  because  the 
solution  of  T  was  with  zero  boundary  conditions;  contributions  from  non-zero 
boundary  conditions  can  be  linearly  superimposed. 

2.  The  method  Is  general  such  that  any  choice  of  basis  functions  can  be 
used.  Also  the  Karhunen-Loeve  condition  Is  imposed.  There  are  also  no 
matrix  Inversions.  Fourier  transform  techniques  deal  with  stationary 
processes  whereas  this  method  results  In  non-statlonary  nonlsotroplc 
covariance  of  the  gravity  field. 
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W.  G.  Heller 

The  Analytic  Sciences  Corporation 
One  Jacob  Hay 
Reading,  MA  01867 


ABSTRACT 

Data  reduction  from  an  airborne  gravity  gradiometer  survey  system  (GGSS)  has 
been  long  recognized  as  a  challenging  problem.  The  difficulty  arises  because  of 
the  large  amount  of  gradient  data  required  to  estimate  gravity  disturbances  at  the 
surface.  The  work  reported  here  decomposes  the  airborne  survey  data  reduction 
problem  into  a  track-by-track  treatment  which  is  both  optimal  in  a  least  squares 
sense  and  utilizes  surface  gravity  vector  tie-point  data.  For  low  altitude  GGSS 
^^urveys  (circa  600m),  it  is  shown  that  the  upward  and  downward  continuation  portions 
of  the  reduction  process  can  be  effectively  accomplished  with  few  data  points 
(e.g.  less  than  10)  and  that  track-by-track  processing  in  a  template  zone 
format  is  near  optimal.  In  this  case,  near-optimal  indicates  that  processing 
algorithm-induced  errors  are  less  than  0.2  mgal/O.OA  arcsec.  Synthetic 
airborne  GGSS  data  from  a  gravity  field  realization  typical  of  the  Clinton-Sherman 
test  area  is  reduced.  Gravity  vector  estimation  results  Indicate  that  surface 
recovery  at  an  accuracy  near  1.0  mgal/0.2  arcsecond  can  be  expected  from  a 
100  E^/Hz  (30  E^/r/sec)  GGSS.  Quantification  of  the  relative  information 
contents  of  each  track  of  gradiometer  data,  as  determined  by  the  gradiometer 
noise,  vs  the  predictability  of  the  gravity  field  along  a  track,  based  on  its 
correlation  characteristics  only  (l.e.  in  the  absence  of  gradient  data), 

'leads  to  further  simplification  in  the  data  reduction  process.  In  particular, 

^j^or  a  gravity  field  typical  of  the  Clinton-Sherman  test  area,  a  statistical 
gravity  model  can  be  dropped  from  the  GGSS  data  reduction  process  idLthout 
penalty  if  the  GGSS  errors  are  less  than  about  300  E^/Hz. 
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TASC'S  ACTIVITIES  IN  THE 
GGSS  DEVELOPMENT  PROGRAM 
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SURFACE  VS  AIRBORNE  GGSS  DATA  REDUCTION 
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ILLUSTRATION  OF  AIRBORNE  GGSS  DATA 

REDUCTION 

OPTIMAL  ESTIMATION  of  T,  (COVARIANCE  ANALYSIS) 
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PAPER  TITLE: 


GRAVITY  GRADIOMETER  (GGSS)  TEST  PLANNING  AND  TEST  DATA  TREATMENT 


SPEAKERS  NAME:  Warren  G.  Heller 
Questions  and  Comnents: 

Stan  Jordan:  You  have  presented  a  simplistic  and  optimistic  view  of  how  tie 
points  are  used.  I  feel  that  the  use  of  tie  points  should  be  approached 
cautiously;  local  geologic  features  may  have  significant  (one  mllligal) 
effects  on  the  tie  points  but  are  unobserved  by  the  gradlometer  at  aircraft 
altitude.  If  such  aliasing  effects  arise  at  the  tie  points,  they  cannot  be 
used  to  Improve  estimation  of  low-frequency  gravity  disturbances  in  the  region. 
These  difficulties  are  not  revealed  by  your  analysis  because  the  7-shell 
AWN  model  is  not  valid  at  wavelengths  shorter  than  3.5  km.  To  avoid  aliasing 
effects  at  tie  points,  mini-surveys  should  be  conducted  in  the  vicinity  of 
each  tie  point.  These  mini-surveys  can  be  done  with  the  (land-based) 
gradlometer  van. 

Julian  Center:  To  do  upward  and  domward  continuation  you  need  to  know  the 
vertical  gradients  all  along  the  path.  This  may  not  be  achievable  with 
high  accuracy  if  there  is  high-frequency  activity  of  the  gravity  field 
around  the  tie  point  or  estimation  point. 

SPEAKERS  RESPONSE:  For  sufficiently  short  upward  continuation  distance 
(e.g.  600  meters),  local  gradient  data  (measured  by  the  gradlometer)  at 
altitude  suffices.  In  an  extreme  case,  say  in  which  the  tie  points  were 
chosen  injudiciously,  such  that  the  gravity  field  curved  significantly 
(perhaps,  SOE  or  more)  then  surface  gradient  measurements  near  the  tie  point 
could  be  appropriate.  The  practicalities  of  continuous  use  of  gradient 
data  along  the  entire  upward  continuation  path  preclude  "continuous  upward 
continuation".  Fortunately  it  is  not  necessary. 
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GRADIENT  INTEGRATION  PROCEDURE  FOR 
PATH  ERROR  REDUCTION 


A.  E.  Rufty 

Naval  Surface  Weapons  Center 
Dahlgren,  Virginia  22448-5000 


ABSTRACT 

A  procedure  to  integrate  gravity  gradients  aloft  to  obtain  gravity 
disturbance  values  (aloft)  will  be  presented.  The  algorithm  simultaneously 
weights  all  possible  path  integrals  of  the  gradient  to  yield  an  optimal 
estimate  of  the  disturbance. 
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OPTIMUM  CHOICE  OF  WEIGHTS 


ASSUMPTIONS 


RESULTS  k^.Smgal) 


PAPER  TITLE:  GRADIENT  INTEGRATION  PROCEDURE  FOR  PATH  ERROR  REDUCTION 

SPEAKERS  NAME:  Alan  E.  Rufty 
'Questions  and  Comments: 

Alan  Zorn:  Stokes'  theorem  provides  a  connection  between  cross  products 
of  vector  fields  and  line  integrals  of  gradients: 

Thus,  the  explicit  constraint  V  x  f^  *  0  (1 e-conservative  field) 

Is  equivalent  to  zero  closed  line  Integrals.  Is  the  explicit  constraint 
really  required  for  your  solution? 

SPEAKERS  RESPONSE:  1  use  It  because  It  Is  easy  to  Implement  In  my  algorithm, 
and  It  does  reflect  the  physics  of  the  problem. 


NASA  REQUIREMENTS  FOR  A  SPACEBORNE  GRAVITY  GRADIOMETER— AN  OVERVIEW 

C.  J.  Finley 

Office  )  Space  Science  and  Applications 
NASA  Headquarters 
Washington,  D-C.  20546 

D.  E.  Smith 

Space  and  Earth  Science  Directorate 
NASA/Goddard  Space  Flight  Center 
Greenbelt,  MD  20771 

ABSTRACT 

Beginning  with  Sputnik  and  Vanguard  In  1958,  NASA  has  used  artifical 
satellites  and  space  platforms  to  measure  and  map  the  Earth's  gravity.  Until 
now,  the  gravity  field  has  been  derived  indirectly  by  the  measurement  of 
satellite  perturbations,  such  as  with  the  Beacon  Explorers  and  Lageous  satellites; 
or  by  means  of  satellite  alitimetry,  such  as  with  GEOS-3  and  Seasat.  Each  of 
these  techniques  suffer  severe  limitations,  the  first  in  resolution  due  to 
the  high  orbital  altitude  of  anywhere  from  500  km  to  almost  6000  km  for  Lageos 
and  the  second  in  both  resolution  and  global  coverage  because  of  the  approximately 
800km  altitude  of  the  altimetrlc  satellites  and  being  effective  only  over 
oceanic  areas.  Spaceborne  gradiometers  will  for  the  first  time  permit  direct 
measurement  of  the  gravity  field  over  the  entire  surface  of  the  Earth  at  much 
Improved  resolutions  and  accuracy.  The  present  gravity  field  with  an  uncertainty 
of  5  to  10  mgals  at  resolutions  greater  than  500  km  is  Inadequate  for  NASA’s 
present  operational  and  scientific  needs.  We  need  1  to  2  mgals  at  100  km 
resolution  for  our  solid  earth  geodynamics  studies  as  well  as  for  Improved 
orbit  determination  and  the  derivation  of  a  10  cm  ocean  geoid  to  be  used  in 


the  Topex  satellite  ocean  circulation  studies. 


In  addition  to  the  Earth,  there  is  equal  or  greater  value  in  being  able  to 
measure  the  gravity  fields  of  the  planets.  The  same  basic  geophysical  questions 
about  the  structural  rigidity  of  the  lithosphere  and  the  existence  and  form  of 
convection  in  the  mantle  apply  to  the  moon  and  the  planets  as  they  do  the  Earth. 
Further,  the  engineering  difficulties  of  carrying  out  a  planetary  mission  compared 
to  an  earth  mission  make  it  essential  that  the  technology  that  is  used  is  fully 
capable  of  meeting  the  scientific  measurement  objectives  in  both  accuracy  and 
resolution;  a  re-flight  may  not  be  a  possibility.  A  space  proven  gravity  gradlometer 
could  be  flovm  on  many  of  the  future  planetary  missions  and  provide  a  wealth  of  new 
information  about  the  planets  and  their  satellites  that  cannot  be  obtained  any  other 
way. 


•  EARTH  STRUCTURE 

•  TECTONIC  MOVEMENT 

•  OCEAN  DYNAMICS 
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NASA  REQUIREMENTS  FOR  A  SPACEBORNE  GRAVITY  GRADIOMETER-AN  OVERVIEW 


SPEAKERS  NAME;  Werner  D.  Kahn 
Questions  and  Comments; 

Jim  Faller:  Given  the  time  slippage  that  has  occurred  with  the  GRM  mission, 
wouldn't  it  be  appropriate  at  this  time  to  once  again  take  a  hard  and  con¬ 
temporary  look  at  the  GRM  mission,  and  in  particular,  to  look  at  what  improve¬ 
ments  could  be  effected  if  a  laser  tracking  system  (perhaps  using  a  stabilized 
laser  diode  light  source)  were  added  and  the  DISCOS  appropriately  upgraded. 

A  laser  tracking  system  could  result  in  two  orders  of  magnitude  tracking 
precision  improvement  which  could,  for  example,  be  used  to  fly  at  a  somewhat 
higher  altitude  (keeping  ground  resolution  the  same)  and  thereby  permit  a 
longer  mission  and/or  a  smaller  vehicle  (which  would  presumably  cost  less). 
Such  an  optimized  mission  should  then  be  compared  to  a  single  satellite 
gradiometer  system.  The  fundamental  difference  between  these  two  approaches 
is  that  one  is  in  effect  a  "point"  Instrument  and  as  such  suffers  terribly 
from  local  internal-to-the-satellite  mass  redistributions  and  their  effects 
in  the  measured  gradlent(s).  The  other,  GRM,  enjoys  the  freedom  from 
these  critical  systematic  error  sources  as  a  result  of  its  long  baseline. 

(End  effects  are  greatly  reduced).  One  should  worry  whether  or  not  the 
single-satellite  gradiometer  can  ever  reduce  the  local  mass  Induced  gradient 
changes  much  below  the  0.01  to  0.1  E  region.  Certainly  to  achieve  10“^  E 
is  difficult,  if  not  impossible,  and  to  achieve  lO'^g  on  the  shuttle  is  simply 
out  of  the  question. 

Since  the  science  to  be  learned  depends  entirely  on  the  achieved  resolu¬ 
tion  on  the  ground  in  the  presence  of  any  and  all  systematic  biases,  it  is 
critical  to  make  realistic  estimates  of  their  magnitudes  and  thereby  the 
useful  as  opposed  to  the  theoretical  sensitivity  of  both  of  these  systems. 

It  is  the  achieved  rather  than  the  theoretical  sensitivity  which  will 
dictate  the  quality  of  the  science.  It  is  my  impression  that  (from  a 
fundamental  design  point  of  view)  an  extended  (GRM)  instrument  (though 
complicated  by  requiring  two  satellites)  enjoys  a  great  advantage  over  the 
simpler  single  satellite  instrument  in  regard  to  this  all  important  and 
pivotal  design  question. 

Klaus-Peter  Schwarz:  The  question  addresses  the  interpretation  of  the 
accuracy  curves.  You  claim  an  accuracy  of  better  than  1  mgal  for  gravity 
anomaly  resolution  and  of  better  than  1  cm  for  geoid  resolution.  On  the 
other  hand,  you  state  a  wavelength  resolution  of  50  km.  These  two  sets  of 
numbers  are  Incompatible  because  the  short  wavelengths  are  not  resolved, 
and  contain  much  large  RMS  errors  than  the  ones  you  quoted. 

SPEAKF.RS  RESPONSE:  The  numbers  refer  to  one  by  one  degree  mean  values.  Short 
wavelength  variations  were  not  considered.  The  RMS-values  do  not  therefore 
represent  gravity  field  or  geoid  resolution  at  individual  points. 
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GRAVITY  FIELD  FINE  STRUCTURE  MAPPING  USING 
A  SPACEBORNE  GRAVITY  GRADIOMETER 

W.  D.  Kahn 

Goddard  Space  Flight  Center 
Mall  Code  621 
Greenbelt,  Maryland  20271 


ABSTRACT 

Covariance  Error  Analysis  studies  have  been  performed  to  assess  the  capabilities 
of  a  Spaceborne  gravity  gradiometer ' s  capability  to  determine  the  short  wavelength 
components  of  the  earth’s  gravity  field.  For  the  studies,  consideration  was  given  to 
the  effects  of  gravity  signal  attenuation  due  to  increase  in  spacecraft  orbit  altitude, 
the  effects  of  orbit  errors,  altitude  rate  errors,  gradiometer  system  precision, 
|estimatlon  strategy  etc.  Typical  results  obtained  are  those  for  a  satellite  in  a  160 
km  polar-circular  orbit,  carrying  a  three-axis  gradiometer  with  a  precision  of  10”3 
EOTVOS  Units  (E)*.  The  satellite  orbit  error  was  assumed  to  be  5.5  meters,  the 
measurement  sampling  rate  was  1  measurement  per  4  seconds.  The  size  of  the  region 
surveyed  over  a  180  day  period  was  19°  x  19°  which  was  traversed  by  approximately  300 
orbit  passes  of  which  a  maximum  of  16  passes  traversed  each  1°  x  1°  gravity  anomaly 
block  within  the  region.  Results  obtained  from  the  simulation  indicate  that  the 
gravity  gradiometer  can  determine  gravity  anomalies  with  a  horizontal  resolution  of 
50  km  to  an  accuracy  of  0.4  mgals.  Slightly  improved  results  say  on  the  order  of  10 
to  20%  can  be  anticipated  if  the  gradiometer  precision  is  Increased  to  4  x  10“^  E. 


*  IE  =  10~^  cm/S/cm 
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•  CIRCULAR/POLAR 

•  ALTITUDE:  200  km 

•  MISSION  DURATION:  180  DAYS 
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•  PRECISION:  10  3E 

•  SAMPLE  RATE:  1  MEAS/4S 

‘  •  SATELLITE  POS.  ERROR:  5.5  METERS 
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PAPER  TITLE: 


GRAVITY  FIELD  FINE  STRUCTURE  MAPPING  USING  A  SPACEBORNE  GRAVITY 
GRADIOMETER 


^^PEAKERS  NAME:  Werner  Kahn 
Questions  and  Cointnents  : 

Alan  Zorn;  Have  you  considered  the  effects  of  local  or  near-field  gradients 
on  your  E  gradlometer ,  especially  In  a  shuttle-based  experiment? 

^  SPEAKERS  RESPONSE:  No  -  but  it  does  need  to  be  considered. 

Dave  Sonnabend:  1.  GRM  should  be  redesigned. 

2.  Gravity  and  other  problems  would  drown  out  shuttle  measurements. 

A  laser  would  improve  GRM  range-rate  measurement,  but  other  error  sources 
would  limit  the  benefit.  Concerns  about  the  shuttle  are  well  founded;  but 
we  don't  intend  to  do  earth  measurements  from  the  shuttle.  In  a  free 
flyer,  we  have  our  work  cut  out  to  justify  that  we  can  live  with  these 
problems . 

Stan  Jordan:  For  the  GRM  and  satellite  gradlometer  missions,  what  is  the 
track  spacing  of  the  orbit  at  the  equator? 

SPEAKERS  RESPONSE:  Speaker  did  not  know  off  hand. 

Lou  Decker:  What  time  frame  is  being  considered  for  the  Space  Shuttle 
Gravity  Gradimeter  Test? 

^^PEAKERS  RESPONSE:  The  initial  flight  test  of  a  spaceborne  Cryogenic  Gravity 
Gradlometer  utilizing  the  space  shuttle  is  projected  for  the  early  1990 's. 
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SUPERCONDUCTING  GRAVITY  GRADIOMETER  ON  THE  SPACE  SHUTTLE 

S .  H .  Morgan 
J.  R.  Parker 
Program  Development 

NASA,  George  C.  Marshall  Space  Flight  Center 
Marshall  Space  Flight  Center,  AL  35812 


ABSTRACT 

NASA  Is  developing  a  highly  precise  Superconducting  Gravity  Gradlometer 
through  the  efforts  of  a  research  group  at  the  University  of  Maryland.  By 
providing  very  precise  global  mapping  of  the  Earth’s  gravity  field,  this 
gradlometer  will  be  an  Important  element  of  NASA's  Earth  Sciences  Program 
of  the  19908.  The  gradlometer  utilizes  superconducting  technology  to  measure 
very  small  differential  gravity  signals  in  the  presence  of  large  accelerations 
and/or  disturbances.  An  engineering  group  at  the  Marshall  Space  Flight 
Center,  in  conjunction  with  a  study  team  composed  of  other  NASA  centers,  the 
U.S.  Army,  the  U.S.  Air  Force,  and  the  University  of  Maryland,  has  begun  an 
analysis  to  define  an  experiment  to  the  flown  aboard  the  Space  Shuttle.  The 
objective  of  the  analysis  Is  to  establish  the  feasibility  and  to  provide  data 
for  future  detailed  design  of  a  Flight  Experiment  System.  This  paper  contains 
results  to-date  of  the  preliminary  design  of  a  Space  Shuttle  Flight  of  the 
Superconducting  Gravity  Gradlometer. 
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DEVELOPMENT  OF  SUPERCONDUCTING  GRAVITY  GRAD IOMETER 
FOR  SPACE  APPLICATIONS 

H.  A.  Chan 
M .  V •  Moody 
H.  J.  Paik 

Dept  of  Physics  and  Astronomy 
University  of  Maryland 
College  Park,  Maryland  -20742 


ABSTRACT 

Principles  and  Design  of  the  Three-Axis  Superconducting  Gravity  Gradiometer 


will  be  discussed.  Status  of  construction  and  test  will  be  reported. 


I 

Development  of  a  Superconducting  Gra 
Gradiomeier  for  Space  applications 


H  A-  Chan  ^  M.V.  Maodtj  ^  an  ef  H-  I-  Paik 
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PAPER  TITLE: 


DEVELOPMENT  OF  SUPERCONDUCTING  GRAVITY  GRADIOMETER  FOR 
SPACE  APPLICATIONS 


P BAKERS  NAME: 


Vol  Moody 


Questions  and  Comments: 

» 

Warren  Heller:  1.  Is  the  common  mode  acceleration  signal  output  from  the 
gradiometer  available? 

2.  Have  you  performed  an  analysis,  relating  the  dimensional  tolerances  of 
the  six-axis  proofmass  to  the  expected  level  of  cross  axis  error  coefficients 
*  expected? 


SPEAKERS  RESPONSE:  1.  Yes,  there  is  a  redundancy  in  the  detection  of 
linear  acceleration  with  the  gradiometer  and  the  accelerometer. 

2.  No,  we  expect  cross  coupling.  The  coefficients  will  be  measured  and 
calibrated  out  using  a  six-axis  shaker. 


PLATFORM  REQUIREMENTS  AND  ERROR  COMPENSATION  FOR  A  SUPERCONDUCTING 

GRAVITY  GRADIOMETER 

H.  J.  Paik 

Dept  of  Physics  and  Astronomy 
University  of  Maryland 
College  Park,  Maryland  20742 

ABSTRACT 

Requirements  of  the  platform  for  10“^  E  gradiometer  will  be  discussed. 

Platform  stabilization  and  error  compensation  by  use  of  a  six-axis  superconducting 
accelerometer  and  a  six-axis  shaker  will  be  discussed. 
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PAPER  TITLE:  PLATFORM  REQUIREMEHTS  AND  ERROR  COMPENSATION  FOR  A  SUPERCONDUCTING 

CRAPITT  GRADIOMETER. 

SPEAKERS  NAME:  Ho  Jung  Palk 

Questions  and  Comnents: 

Alan  Zorn:  You  promised  earlier  to  discuss  the  effect  of  self-gradients  on 
a  satellite-based  cryogenlc-gradlometer  with  10“^  accuracy.  Is  this 

a  good  time? 

SPEAKERS  RESPONSE:  Yes,  self-gradients  contribute  a  DC  effect,  and  we  do  not 
consider  them. 


Alan  Zorn:  But,  self-gradients  are  not  constant  (eg-out gas sing,  transfer  of 
mass  %rlthln  helium  cooling  system,  etc)  and  would  certainly  appear  at  10~^ 

E/(Hz)1/2 


SPEAKERS  RESPONSE:  He  can  locate  such  moving  items  far  from  the  sensing 
elements.  He  have  been  assured  that  the  helium  cooling  system  will  not 
affect  such  accuracy.  The  next  speaker  may  be  willing  to  further  expound  on 
this  question. 

Warren  Heller:  Please  comment  on  the  practicality  of  your  requirement  to 
control  Jitter  to  10“^^  rad  (or  even  10^®  -  less  stringent  case). 


SPEAKERS  RESPONSE:  With  the  co-rotating  scheme  that  I  proposed,  the  angular 
stability  of  the  dewar  should  be  very  good.  Our  angular  accelerometer  can 
be  tuned  to  measure  angular  acceleration  down  to  10“^^  rad/s^/CHz)^/^.  po^ 
10”^  E/Hz^/2^  required  angular  velocity  sensitivity  is  4x10"^^  rad/s/Hz^/^ 
for  the  worst  case  which  is  for  a  local  geographic  orientation.  This  will  be 
no  problem  for  our  angular  accelerometer. 


TLC  FOR  A  MAGNETICALLY  FLOATED  GRAVITY  GRADIOMETER 


D .  Sonnabend 

Jet  Propulsion  Laboratory 
California  Institute  of  Technology 
4800  Oak  Grove  Drive 
Pasadena,  California  91109 


ABSTRACT 

The  proposed  NASA  Gravity  Gradiometry  Mission,  a  very  low  altitude  satellite, 
carrying  a  cryogenic  gradiometer  under  development  at  the  University  of  Maryland, 
poses  some  unique  technical  challenges.  Like  some  other  highly  sensitive  inertial 
instruments,  this  one  can  deliver  its  promised  performance  only  in  orbit,  and  thus 
cannot  be  fully  tested  in  the  laboratory.  Even  in  orbit  we  face  problems  from  vibra¬ 
tion,  acceleration,  rotation,  and  local  gravity  noise.  Our  approach  to  this  is  based 
on  a  magnetic  isolation  system,  imder  development  at  JPL,  and  in  cooperation  with 
other  efforts  at  Marshall  and  Goddard.  The  paper  will  discuss  how  we  will  float  the 
whole  experimental  dewar,  and  measure  the  Instrument's  altitude  and  operating 
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COMMON  MODE  BALANCING  GRADIOMETER  WITH  MONOCRYSTALLINE  SILICON  SUSPENSION 
FOR  HIGH  SENSITIVITY  GRAVITY  MEASUREMENTS 
Jean-Paul  Richard 

Department  of  Physics  and  Astronomy 
University  of  Maryland 
College  Park,  Maryland  20742 


ABSTRACT 

Noise  In  a  room  temperature  gradlometer  with  a  new  low  frequency  common  mode 
balancing  monocrystal line  silicon  suspension  and  laser  instrumentation  is  analyzed. 
The  potential  sensitivity  as  limited  by  the  noise  associated  with  damping  in  silicon 
and  interferometric  position  sensing  Is  projected  at  2x10”^  Eotvos/  Hz^/2. 
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PAnSR  TITLE: 


OOHHM  MODE  BALARCING  GRADIOMBTBR  WITH  MONOCRTSTALLINE 
SILICON  SUSPENSION  FOR  HIGH  SSNSITIVITT  ^VITT  MEASUREMENTS 


SPEAKERS  NAME:  Jean-Paul  Richard 


esclons  and  Conaents: 


Dave  Sonnabend:  Don't  use  Invar  for  your  structure,  magnetos trlctive  strain 
**  Is  worse  than  temperature  strain  In  routine  materials* 


Andrew  Lazarewlcz:  With  such  a  high  Q,  isn't  your  bandwidth  very  narrow? 


•  SPEAKERS  RESPONSE:  The  mechanical  Q  Is  sharp,  so  is  the  mechanical  bandwidth 
narrow.  The  Instrument  bandwidth  is  not  limited  by  the  mechanical  bandwidth. 
Note  that  both  signal  and  noise  Is  amplified  within  the  mechanical  bandwidth, 
so  S/N  may  not  be  Improved  (notable  exception:  shot  noise).  In  fact,  most 
gravitational  signals  of  Interest  are  outside  the  mechanical  bandwidth  in 
question. 


Jim  Pallet:  Why  not  measure  the  length  of  the  bar  with  an  Independent  laser 
&  thereby  take  out  (servo  out)  temp-effects:  Servo  to  keep  fixed  &  then 
measure  distance  with  a  laser  to  "measure"  gradients. 


LIQUID-SUPPORTED  TORSION  BALANCE  AS  GRADIOMETER 
J.  E.  Faller 
P.  T.  Keyser 

Joint  Institute  For  Laboratory  Astrophysics 
University  of  Colorado 
Boulder,  Colorado  80309 

ABSTRACT 

We  employ  a  liquid-supported  torsion  balance  as  a  fixed-site  long-term  curvature 
variometer.  The  traditional  fiber  is  replaced  by  liquid  support  and  electrostatic 
positioning.  Thus  the  torsion  constant  is  adjustable  by  varying  the  voltage  applied 
to  the  torque  electrodes,  while  the  centering  voltage  remains  constant.  The 
sensitivity  of  this  type  of  gradlometer  will  be  discussed,  along  with  critical 
parameters  for  success.  Preliminary  data  will  be  presented. 


Liquid-Supported  Torsion  Balance  as  Gradloneter 


P.  T.  Keyser  and  J.  E.  Faller 
Joint  Institute  for  Laboratory  Astrophysics, 

University  of  Colorado  and  National  Bureau  of  Standards, 

Boulder,  Colorado,  80309-0440 

We  employ  a  liquid-supported  torsion  balance  as  a  fixed-site  long¬ 
term  curvature  variometer.  The  traditional  torsion  fiber  Is  replaced  by 
liquid  support  and  electrostatic  positioning.  Thus  the  torsion  constant 
is  adjustable  by  varying  the  voltage  applied  to  the  torque  electrodes, 
while  the  centering  voltage  remains  constant.  The  sensitivity  of  this 
type  of  gradlometer  will  be  discussed,  along  with  critical  parameters 
for  success.  Preliminary  data  will  be  presented. 


The  world's  first  gradiometers  were  the  fiber-supported  torsion 
balances  (FSTB's)  used  by  Baron  Roland  von  EOtvOs  in  the  late  nineteenth 
century.^  A  wide  variety  of  supported  mass  configurations  were  devel¬ 
oped  (the  two  basic  ones  are  shown  in  Figure  1)  and  such  torsion  balance 
gradiometers  were  used  In  geophysical  research  through  about  1940,  at 
which  time  gravimeters  replaced  them.  Torsion  balances  themselves  were 
Invented  in  the  late  eighteenth  century  and  have  been,  and  still  are, 
used  in  a  variety  of  applications  for  measuring  very  small  forces.  In 
our  own  work  on  the  Einstein  Equivalence  Principle  we  have  constructed. 


(Gradlooccer) 


(Curvature  Varloaeter) 


FIGURE  1. 

Baron  Roland  von  BStvda*  FSTB  Gradloneters  (1896) 
(Ref.  1,  pp.  369-3701 


ostensibly  as  an  auxiliary  Instrument,  a  pair  of  torsion  balance  gradlo- 
aeters  to  monitor  longterm  (T  ■>  1  day)  changes  In  what  may  be  loosely 
called  the  "horizontal  gravity  gradients"  (more  precisely:  In  the  curva¬ 
tures  of  the  gravity  level  surface).  However,  due  to  the  often  bizarre 
and  always  subtle  problems  associated  with  the  use  of  fibers,  and  due 
to  our  desire  to  Increase  the  sensitivity  of  the  torsion  balance,  we 
have  used  a  radically  different  design.  FSTB's  have  an  Inherently  limi¬ 
ted  sensitivity  to  mass-dependent  forces  as  may  be  seen  by  noting  that 
the  supported  weight  (test  masses  and  beam)  Is  proportional  to  the 
cross-sectional  area  of  the  fiber,  while  the  observed  angle  Is  propor¬ 
tional  to  the  square  of  the  area.  This  results  In  a  maximum  sensitivity 
as  determined  by  the  fiber  material  (yield  strength  versus  elasticity) 
and  by  the  resolution  of  the  detection  system. 

To  overcome  this  sensitivity  limit,  the  supported  object  can  be 
Immersed  In  a  liquid,  to  decrease  the  load  on  the  fiber,  as  John  Henry 
Poyntlng  first  suggested.*^  However,  to  completely  avoid  the  use  of 
fibers,  the  recently-perfected  liquid-supported  torsion  balance  (LSTB) 
uses  an  electrode  array  to  provide  the  centering  and  restoring  torque 
usually  provided  by  the  fiber.®  Our  LSTB  curvature  variometer  is  pic¬ 
tured  In  Figure  2.  An  advantage  of  this  "electrostatic  fiber"  Is  that 
the  voltage  applied  to  the  center  electrode  may  be  (and  usually  Is)  much 
larger  than  that  applied  to  the  torque  electrodes,  which  allows  the  cen¬ 
tering  and  torqulng  forces  to  be  effectively  Independent.  (The  force 
between  two  electrodes  Is  proportional  to  the  square  of  the  applied 
voltage,  and  for  spherical  electrodes  the  force  Is  almost  In  inverse 
proportion  to  the  gap  between  the  upper  and  lower  electrodes.)  Further¬ 
more,  the  torque  voltage  Is  readily  adjustable  so  that  we  can,  for 
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example,  let  the  system  thermally  equilibrate  with  high  torque  voltage 
(“stiff  fiber")  and  then  simply  lower  the  torsion  constant  to  the  cor¬ 
rect  value  for  critical  damping  (the  optimum  operating  point).  Similar¬ 
ly,  we  can  quickly  establish  the  true  “zero"  of  the  system  so  as  to  be 
able  to  measure  the  dc  curvature  values.  A  polyatomic  Inert  gas  (such 
as  well-filtered  N2)  used  as  the  fill  gas  maximizes  the  dielectric 
strength  of  the  electrode  gaps  and  thus  yields  the  maximum  stiffness  of 
the  electrostatic  fiber. 

In  a  LSTB  there  are  only  a  few  undesirable  mechanical  oscillation 
modes,  most  of  which  have  frequencies  much  higher  than  the  torsional 
oscillation  frequency  (T  ••  10  min),  and  all  of  which  are  rapidly  damped 
by  the  liquid.  By  contrast,  FSTB's  are  plagued  with  a  great  number  of 
hlgh-Q  modes  which  readily  couple  to  the  torsion  mode.^  The  bobbing 
mode  of  a  LSTB  (T  "  1  sec)  has  never  been  observed  to  couple,  while  the 
off-centering  mode  (T  20  min)  Is  rarely  excited  and  Is  highly  damped. 
There  are  two  tilting  modes  (T  ■  1  sec)  one  of  which  (or  both.  If  they 
are  degenerate)  will  of  necessity  couple  to  the  torsion  mode.  This  may 
be  most  readily  seen  by  noting  that  a  book  (held  shut  with  a  rubber 
band)  may  be  flipped  Into  the  air  about  any  of  three  principal  axes  - 
and  about  two  of  those  the  rotation  Is  stable.  About  the  axis  with  the 
Intermediate  moment  of  Inertia  (or  about  the  two  axes  of  equal  moments 
of  Inertia)  the  motion  Is  unstable.^  Due  to  various  practical  con¬ 
straints,  In  any  real  LSTB  the  moment  of  Inertia  about  the  torsional 
axis  will  always  be  smaller  than  the  other  two  moments,  and  hence  the 
torsional  oscillation  Is  always  stable.  Moreover,  the  coupling  of  the 
tilting  mode(s)  to  the  torsional  mode  only  results  In  short  episodes  of 
small  high-frequency  "noise"  (and  the  tilt  modes  are  rarely  excited). 


We  use  a  transmission-optics  variant  of  the  Gauss  optical  lever  as 
our  angle  detection  system.  The  Infrared  LED  (noted  In  Figure  3)  is  a 
point  source  and  the  lens  Is  a  (boroslllcate  galss)  rod  %d.th  its  axis 
vertical.  The  rays  from  the  LED  pass  through  the  rod  perpendicular  to 
Its  axis  and  are  focused  Into  a  vertical  line.  The  focal  plane  (l.e., 
the  line  Image)  Is  made  to  coincide  with  the  plane  of  the  split  photo¬ 
diode.  The  two  currents  generated  In  the  photodiode  halves  are  convert¬ 
ed  to  voltages  In  low-current-nolse  preamplifiers  and  their  sum  and  dif¬ 
ference  formed.  The  four  sums  and  four  differences  are  combined  In  ap¬ 
propriate  ways  to  generate  radial  translation  (off-centering)  and  angu¬ 
lar  azimuth  signals.  The  net  gain  Is  such  that  we  have  a  measured  sen¬ 
sitivity  of  about  30  mV/arcsec,  using  4.5  mA  In  the  LED. 

The  second  derivatives  of  U,  the  scalar  gravity  potential,  form  a 
tensor,  whose  nine  components  are  reduced  by  the  constraints  Imposed  by 
the  curl  of  g  and  the  divergence  of  g  to  five  linearly  independent  com¬ 
ponents.  Of  these  five  we  happen  to  be  Interested  only  In  those  two 
which  correspond  to  masses  near  the  Instrument  horizon,  that  la  the  two 
curvature  values.  In  fact  the  azimuthal  torque  on  a  LSTB  Is  given  by 

(1)  N  -  1/2  sin(2a)(Ix*x’  "  Iy’y')(Uyy  -  Uxx) 

where  a  Is  the  azimuth  of  the  Instrument  relative  to  the  axes  of  the 
principal  radii  of  curvature,  Uxx  ^yy  principal  cur¬ 
vature  values  (so,  e.g.,  Rx  “  Ez/^xx)*  the  I's  are  the  two 

horizontal  moments  of  Inertia.  Because  of  the  symmetry  of  our  float 
this  reduces  to  (along  the  sensitive  axes  which  are  at  ±45*  to  the  line 
Joining  the  test  masses): 


(2) 


(■r^XUyy  -  Ujx) 


where  m  Is  the  mass  of  each  test  mass  (2.935  kg)  and  r  is  its  radius 
(8.65  cm).  For  an  external  (point  or  spherical)  mass  M  located  on  the 
sensitive  axis  (x  or  y)  at  a  distance  R  from  the  axis  of  the  torsion 
balance,  we  have  for  the  torque 

(3)  N  -  (mr^)(2GM/R^) 

For  our  apparatus,  with  an  external  mass  M  «  150  kg  at  R  *  2  m  (corres¬ 
ponding  toO*2.5EU,  lEO*  10“^  sec~^)  we  find  a  torque  of  about  0.55 
mdyne-cm.  Our  critically-damped  torsion  constant  is  1.09  dyne-cm/rad 
(at  a  voltage  of  about  55  V  rms,  electrode  gap  about  0.1  in.),  so  that 
the  resultant  angle  is  0.50  mrad  or  104  arcsec.  This  gives  a  sensitivi¬ 
ty  of  over  40.  arcsec/EU,  or  to  put  it  another  way  the  2.5  EO  signal  is 
over  40  times  the  rms  noise  (from  all  sources)  in  our  apparatus. 

The  noise  sources  may  be  conveniently  characterized  as  Internal  or 
uncorrelated  (l.e.,  imcorrelated  with  any  external  signal)  and  external 
or  correlated.  Tests  have  shown  that  the  external  noise  sources  have  a 
very  small  effect. 

The  noise  is  observed  to  be  insensitive  to  atmospheric  pressure 
(the  tank  is  sealed  and  fairly  rigid).  We  have  carefully  avoided  mater¬ 
ials  (or  inclusions)  of  high  magnetic  susceptibility,  and  we  keep  all 
large  magnets  well  away  from  the  apparatus,  so  that  the  only  significant 
noise  of  magnetic  origin  is  due  to  eddy  currents  produced  by  rotating 
magnetic  fields.  It  is  quite  possible  to  produce  rotating  fields  of 
sufficient  frequency  and  Intensity  to  cause  the  float  to  rotate  (at  up 


to  1  rpm),  and  we  do  so  on  occasion,  but  the  ordinary  60  Hz  rotating 
magnetic  fields  produced  by  lab  motors  have  no  observable  effect.  (In 
addition,  we  have  obtained  62-ratl  high-permeabllity  MIL-N14411C  80%-Ni 
magnetic  shields  which  we  plan  to  Install.)  Ordinary  floor  vibrations 
(people  jumping  at  R  ■  2  m,  e.g.)  have  no  observable  effect,  but  large 
earthquakes  (as  in  Figure  4)  or  pathological  (badly  unbalanced)  machines 
(as  in  Figure  5)  which  produce  large  horizontal  accelerations  or  large 
tilts  of  the  floor  do  have  an  effect.  Such  noise  is  intermittent  and 
readily  identifiable,  and  causes  no  real  difficulty. 

Internal  noise  sources  cause  more  difficulties,  none  Insurmount¬ 
able.  The  torque  and  center  voltages  are  ac  to  avoid  the  possibility  of 
.Jtatlc  charge  accumulation  (though  all  exposed  surfaces  of  the  float  are 
conducting  and  grounded,  and  tests  with  dc  voltages  never  showed  any  ef¬ 
fects  attributable  to  static  charge).  The  center  voltage  Is  regulated 
to  AV/V  ■  10“^  and  the  torque  voltage  to  AV/V  ■  SxlO"**.  The  (spurious) 
change  of  angle,  A9,  that  arises  from  a  shift  In  torque  voltage  Is 
1O“^0.  The  observed  (gausslan)  voltage  noise  In  the  optical  detection 
system  amounts  to  15  mV  (pp)  -  l.e.,  0.5  arcsec. 

The  dominant  noise  sources  are  due  to  convection  currents  In  the 
liquid,  which  may  arise  either  from  Impurity  concentration  gradients  or 
from  thermal  gradients.  Nothing  Is  more  crucial  to  the  success  of  a 
LSTB  than  the  elimination  of  these  two  problems. 

Water  Is  the  universal  solvent,  and  we  go  to  great  lengths  to  pur¬ 
ify  the  water.  We  use:  predelonizatlon,  followed  by  5um  filtration,  UV 
sterilization,  I  ym  filtration,  carbon  filtration,  deionization,  and 
0.22  ym  filtration.  All  but  the  first  stage  are  In  a  recirculating 


system 


6  Sept  1985,  03:47:29  GMT 
43.2“N,  110. 8"W;  R  -  4.6 
(near  Jackson,  Wyoming) 


FIGURE  4. 

Effects  of  Earthquakes  on  LSTB 


400 
uVol  t 


HftCHntE  Cit<  ■H^S> 


FIGURE  5. 


Horizontal  Floor  Accelerations  (2  uV/ym) 


The  resistivity  of  acceptable  water  is  10  MO-ca  (though  in  practice  we 
usually  obtain  15  MO-cn).  Experience  shows  that  it  is  necessary  to  add 
a  liquid  non-ionic  surfactant  to  the  water  to  form  a  liquid  surface  mon¬ 
olayer  which  displaces  the  otherwise  inevitable  solid  surface  layer 
formed  of  residual  (less  than  1  ppm)  impurities.  We  use  Kodak  Photoflo 
2100  (not  the  premix  Photoflo  200)  or  Triton  KlOO,  at  about  500  ppm.  In 
addition  it  is  crucial  to  protect  the  aluminum  from  corrosive  attack  by 
the  water.  (Such  attack  is  aggravated  by  dissolved  oxygen,  so  we  deoxy¬ 
genate  the  water  by  bubbling  filtered  N2  through  it.)  The  main  defense 
that  we  use  against  this  otherwise  Inevitable  corrosion  is  the  applica¬ 
tion  of  a  commercial  chromate-conversion  r mating  to  the  well-cleaned  and 
unetched  aluminum  (as  per  MIL-C554l).^  We  do  not  dissolve  chromate  ions 
in  the  water  (which  has  also  been  shown  to  prevent  corrosion  effective¬ 
ly)  as  there  is  some  indication  that  they  react  with  our  surfactants. 
Without  this  protective  coating  the  LSTB  falls  (“locks  up"  in  a  weak 
surface  gel)  in  about  a  week,  but  with  the  coating  we  can  operate  for 
eight  months  or  more. 

Water  is  used  because  it  is  cheap,  safe,  readily  purified,  and  be¬ 
cause  it  alone  of  all  liquids  has  a  maximum  density  point. ^  Most  therm¬ 
al  convection  currents  in  a  gravity  field  are  driven  by  buoyancy  and 
these  bouyancy  forces  arise  due  to  the  temperature-dependent  density  of 
the  liquid  (similar  to  the  principle  of  the  hot-air  balloon).  At  3.98 
®C  the  expansion  coefficient  of  water  is  zero  —  so  that  small  thermal 
gradients  do  not  cause  convection  currents.  There  are  five  distinct 
sources  of  thermal  gradient  forces  on  a  float.  Two  of  these  forces 
(which  arise  from  the  temperature-dependance  of  the  surface-tension  and 
the  pressure)  decenter  the  float  and  need  not  concern  us  here.  The 


other  three  forces  (which  are  due  to  various  types  of  convection  cur¬ 
rents)  can  torque  the  float.  The  dominant  torquing  convection  current 
is  the  circumferential  current  in  the  water  between  the  vertical  walls 
of  the  float  and  tank.  This  is  analogous  to  the  "double-pane  window" 
problem,  where  the  idea  is  to  ralniraljte  the  heat  flux.^®  Here  however  we 
seek  to  minimize  the  convection  (it  can  never  go  to  zero  unless  the 
temperature  gradient  or  the  width  of  the  gap  is  exactly  zero).  In  the 
case  of  this  circumferential  current  the  maximum  allowable  AT  (across 
the  tank)  for  a  specified  Ad  (at  critical  damping)  scales  as: 

(4)  AT  «  (r2  -  ri)“^ 

so  that  (all  else  being  equal)  we  should  minimize  this  gap. 

We  go  to  great  lengths  to  insulate  and  control  the  temperature  of 
our  apparatus.  A  nultilayer  cubical  enclosure  (alternating  insulation 
and  reflective  metal  cans)  three  feet  on  a  side  surrounds  the  one-foot 
cylindrical  tank.  Around  the  innermost  of  these  outer  cans  is  wrapped 
flexible  hose  through  which  temperature-controlled  water  flows.  Around 
this  same  can  is  wrapped  a  heater  connected  to  a  proportional-integral 
controller  (gain  120  dB,  time  constant  11000  sec,  stability  better  than 
0.1  m®C,  resettablllty  1  m®C).^^  The  temperature  of  this  composite  sys¬ 
tem  is  set  to  give  4.0  ±0.1  “C  at  the  waterline  of  the  tank.  (A  seren¬ 
dipitous  effecc  of  operation  at  4  ®C  is  that  the  corrosive  attack  of 
water  is  slowed  by  over  two  orders  of  magnitude  from  its  rate  at  20 
®C.)  The  net  effect  is  a  stability  of  better  than  0.5  m°C  on  the  tank 
over  several  days  (Figure  6)  and  a  sensitivity  to  room  temperature  fluc¬ 
tuations  of  less  than  2  m®C/'’C.  We  believe,  in  spite  of  all  this,  that 
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FIGURE  6. 


Tank  Temperature  at  the  Waterline 
(Start;  850915,  22:13) 


our  dominant  noise  term  Is  temperature-induced  convection  currents. 

In  Figure  7  we  present  sample  data  from  a  run  with  an  Identical 
system  In  which  the  two  masses  were  replaced  %flth  a  disk  (of  the  same 
total  mass  and  moment  of  Inertia  about  z).  In  such  a  device  no  torque 
should  arise  from  curvatures  and  so  any  signal  Is  "noise'*  -  torques  due 
to  floor  vibrations  (the  10-mln  period)  or  torques  due  to  convection 
currents.  Note  that  In  two  weeks  (16  Jan  to  1  Feb  1986)  the  "zero"  has 

drifted  almost  0.2  V  or  0.4  arcsec/day  or  0.01  EU/day.  Some  sample  data 

from  the  gradlometer  system  are  presented  In  Figures  8  through  11.  As 
can  be  seen,  a  claimed  sensitivity  of  0.1  BU  would  not  be  amiss.  The 
drift  rate  In  Figure  11  Is  less  than  75  mV/day  or  2.5  arcsec/day  or 

0.062  EU/day.  (This  drift  was  measured  during  a  period  we  now  have  rea¬ 

son  to  suspect  marked  the  beginning  of  the  degradation  of  the  water,  af¬ 
ter  eight  months  of  operation.  Before  this  It  seemed  to  be  about  5 
times  less,  consistent  with  the  data  of  Figure  7.)  It  must  also  be 
noted  that  the  response  time  Is  vastly  slower  than  the  various  dynamic 
gradlometers  now  available.  In  fact  the  damping  time  Is  almrst  20  min 
(so  that  the  response  to  a  step  will  be  95Z  complete  after  about  three 
such  times  or  one  hour)  which  Is  about  two  orders  of  magnitude  slower. 

Future  plans  Include  the  construction  of  three  new  machined  float- 
and-tank  sets  (the  current  sets  are  modified  spun-aluminum  cook  pots). 
These  new  sets  will  have  a  0.5  cm  gap  between  the  float  and  tank  which 
should  reduce  the  effect  of  the  convection  currents  by  a  factor  of  about 
100,  and  thus  floor  vibrations  and  electronics,  not  convection,  should 
be  the  dominant  noise  sources.  In  addition,  the  larger  tungsten-alloy 
masses  will  be  at  a  greater  radius  for  a  gain  of  a  factor  of  2.64  In 
sensitivity  to  curvatures.  The  first  set  Is  ready  and  will  be  given  a 


preliminary  test  as  soon  as  the  larger  masses  are  ready. 

This  work  was  supported  by  the  National  Bureau  of  Standards,  the 
Sensor-Technology  Division  of  the  Belvolr  Research  Development  and  Engi¬ 
neering  Center,  and  the  Air  Force  Geophysics  Laboratory. 
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FIGURE  9. 

S«all  Signal  (8  Sept  1985) 
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FIGITRE  11. 

Float  Position  and  Fourier  Transform 
(Start  850819,  18:16) 
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LIQUID-SUPPORTED  TORSION  BALANCE  AS  GRADIOMETER 


SPEAKERS  NAME:  Paul  T.  Keyser 
Questions  and  Comments : 

Peter  Ugincius:  Would  your  instrument  be  able  to  detect  solid-earth  tides? 

SPEAKERS  RESPONSE:  Don't  know.  Depends  on  what  U^^,  are  for  solid-ea 

tides. 
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ABSTRACT 

The  development  of  a  one-axis,  non-cryogenic  gravity  gradiometer  is 
described.  The  Instrument  presently  developed  at  IFSI/CNR  (Frascati,  Italy) 
consists  of  two  displacement  sensors,  50  cm  apart  in  the  instrument  prototype. 
Since  the  sensitive  axes  of  the  sensors  are  orthogonal  to  the  baseline  the 
gradiometer  measures  the  offline  components  of  the  gravity  tensor.  The  proof 
^^sses,  rectangular  in  shape,  are  connected  to  the  main  frame  by  a  crank 
shaped  attachment;  every  part  being  machined  integrally  from  the  same  A1 
plate  which  constitutes  the  rigid  baseline  of  the  instrument.  Displacements 
of  each  proof  mass  are  measured  by  capacitive  transducers  whose  outputs  are 
sent  to  an  FET  amplifier  with  noise  temperature  Tn  =  100  mk.  A  feedback  loop 
provides  a  constant  electric  field  on  the  moving  plates,  resulting  in  a 
dynamic  range  enhancement  and  sensitivity  increase.  A  sensitivity  of  the 
order  of  10~2  EU/Hz^/^  ig  expected,  sufficient  to  perform  significant  measurements 
of  the  eartn’s  gravity  field  from  a  low  orbit  spacecraft  such  as  TSS-2. 

Results  from  preliminary  laboratory  tests  of  the  Instrument  are  also  illustrated. 
'^In  the  ground,  more  exhaustive  tests  will  be  performed  jointly  with  the 
Smithsonian  Astrophyslcal  Observatory  by  using  a  free  fall  technique  in 
\accuum  in  order  to  isolate  the  Instrument  from  seismic  and  acoustic  noise. 
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RESULTS  OF  LABORATORY  TESTS  PERFORMED  SO  FAR  INDICATE  THAT  THE  PREDICTED 
INSTRUMENT  PARAMETERS  ARE  ACHIEVABLE  IN  PRACTICE.  THEREFORE  THE  EXPECTED 


PAPER  mLE: 


DEVELOPMENT  OF  A  HIGH  SENSITTVITf,  NON-CRYOGENIC  GRAVITY  GRAPTOMETER 


tPEARERS  NAME:  Enrico  FuUgnl 
Questions  and  Comments: 

Dave  Sonnabend :  I  don't  believe  nuisance  dynamic  angular  rate  in  falling 
“  elevator  can  be  measured  to  adequate  accuracy. 

|SPEAKERS  RESPONSE:  Not  sure. 

^Warren  Heller:  Common  mode  acceleration  must  be  measured  and  controlled  to 
one  part  in  ten  to  the  ninth.  I  did  not  see  that  this  iss  le  was  nddressed. 

Ho  Jung  Paik:  1.  Your  basic  design  involving  only  two  proof  masses  makes 
your  device  an  angular  accelerometer.  At  10~^E  level,  you  need  to  compensate 
an  angular  acceleration  to  10“^^  rad  s~^.  Which  gyroscope  would  you  use  to 
compensate  this  noise? 

2.  Even  if  there  is  no  angular  acceleration,  there  will  be  a  centrifugal 
acceleration  due  to  angular  velocity. 

SPEAKERS  RESPONSE:  1.  There  is  no  angular  acceleration  because  the  instrument 
will  be  tested  in  free  fall. 

2.  The  geometry  of  the  instrument  makes  it  insensitive  to  centrifugal 
acc<*leration . 
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Ho  Jung  Paik:  That's  not  true.  If  you  have  an  angular  velocity  45°  away 
rom  the  baseline,  you  will  be  sensitive  to  the  centrifugal  acceleration, 
is  will  be  true  even  if  you  employ  4  masses  as  Metzger  does.  It  has  been 
proven  that  the  centrifugal  acceleration  cannot  be  taken  out,  even  in  principle 
in  a  second-order  gradiometer,  without  stabilizing  the  platform. 
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